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CHAPTER  1 
INTRODUCTION 


Nuclear  magnetism  in  metals  has  been  a  very  rich  field  of  study  for  experimental 
and  theoretical  physics.  These  systems,  however,  cannot  be  studied  using  standard 
laboratory  apparatus.  The  small  size  of  the  magnetic  moments,  and  their  weak 
interactions,  call  for  experiments  at  ultra-low  temperatures  with  extremely  sensitive 
probes. 

The  characteristics  of  the  noble  metals  (Cu,  Ag,  Au)  makes  them  ideal  candidates 
for  studying  nuclear  magnetism  in  metals.  The  interactions  of  the  nuclear  moments 
in  these  systems  are  well  characterized,  and  experiments  [1,  2,  3]  on  these  systems 
provide  a  clean  test  of  many  interesting  magnetic  phenomena. 

As  compared  to  electronic  moment  materials,  the  interactions  in  these  systems 
are  very  small.  By  scaling  the  nuclear  magneton  to  the  Bohr  magneton  (which  char- 
acterizes electronic  magnetism),  one  finds 

=  I2»  „  iQ-3,  (1.1) 

so  that  the  interactions  between  magnetic  moments  (which  vary  as  fi2)  are  about  six 
orders  of  magnitude  smaller  than  for  electronic  momemts.  While  the  temperature 
range  of  interest  for  electronic  magnets  is  typically  on  the  order  of  several  Kelvins 
and  higher,  for  nuclear  magnets  one  has  to  work  at  fiK  or  lower  temperatures.  This 
leads  to  the  possibility  of  studying  the  nuclear  spins  not  only  in  their  ordered  state 
[2],  but  also  at  negative  spin  temperatures  [3],  due  to  their  long  spin-lattice  relaxation 
times. 
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1.1    Nuclear  Magnetism  in  the  Noble  Metals 


The  Hamiltonian  H  for  the  nuclear  spins  in  the  noble  metals  is  particularly  simple, 
and  enables  comparison  between  theory  and  experiments.  In  finite  magnetic  field  B0, 
it  consists  of  the  dipolar,  indirect  exchange  (RK)  [5],  and  Zeeman  interactions  between 
the  nuclear  moments.  The  three  terms  are  of  the  form 

HdipolaT  =  ^E(^[£  £  -  3(4  ■  fy)C5  •  Ml)  ,  (1-2) 

47r    i<j  \  rij  ' 


i<j 


iqif(2kFrij)ii  ■  ij 


m3 
ij 


(1.3) 


and 


Hz  =  -hB0-J2^u  (1-4) 


where  /x0  is  the  permeability  of  free  space,  h  is  the  Planck  constant,  7*  is  the  gyromag- 
netic  ratio  of  the  nuclear  isotope,     is  the  distance  between  the  spins,  kp  is  the  Fermi 

wave-number,  and  the  /j's  are  the  nuclear  spin  operators.  The  relative  magnitude  of 
the  RK  and  dipolar  interactions  is  characterized  by  the  parameter  77.  As  the  noble 
metals  all  have  a  face-centered-cubic  structure,  one  can  take  the  Fermi  surface  to  be 
spherical.  This  approximation  gives 

f(2kFrij)  =  cos(2kFrij)  -  sm(2kFrij)/(2kFrij).  (1.5) 

for  the  distance  dependence  of  the  interaction  between  nuclear  spins.  The  oscillatory 
behavior  of  the  RK  interaction  can  lead  to  ferromagnetic  or  antiferromagnetic  inter- 
actions between  the  nuclear  spins.  The  cubic  symmetry  of  the  lattice  excludes  any 
electric  quadrupolar  (or  higher  order)  interactions  between  the  moments  and  electric 
field  gradients  in  the  system  [6]. 
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Table  1.1:  Nuclear  magnetic  properties  of  the  noble  metals. 


Isotope 

7/2?r  (kHz/mT) 

k  (K-s) 

T-2  (ms) 

197  Au 

0.143 

0.741 

4.6 

2 

107  Ag 

-0.113 

1.812 

12 

10 

109  Ag 

-0.130 

2.084 

9 

10 

63Cu 

2.22 

11.55 

1.27 

0.15 

65Cu 

2.38 

12.38 

1.09 

0.15 

In  Table  1.1  are  listed  several  of  the  relevant  physical  properties  [7]  of  the  noble 
metals:  the  magnet  moment  fj,,  the  gyromagnetic  ratio  7/(271-),  the  Korringa  constant 
k,  and  the  spin-spin  relaxation  time  T^. 

1.2    Nuclear  Magnetic  Ordering  in  Metals 

When  the  nuclear  spin  temperature  is  on  the  order  of  the  interaction  energy  of 
the  spins,  it  is  possible  for  the  moments  to  become  magnetically  ordered.  The  phase 
diagram  in  copper  [2]  and  silver  [8,  9]  has  been  studied  at  positive  temperatures,  where 
the  competition  between  the  anisotropic  dipolar  interaction  and  the  RK  interaction 
lead  to  a  rich  set  of  ordered  structures  [10].  The  transition  temperature  for  copper  is 
found  to  be  56  nK  [11],  while  for  silver  it  is  560  pK  [12].  Experimentally,  it  is  found 
that  silver  is  an  antiferromagnetic  exchange  dominated  system  [3],  while  in  copper 
the  nuclear  spin  interaction  is  dominated  by  the  anisotropic  dipolar  interaction  [1]. 

At  constant  postive  temperatures,  the  ordered  arrangement  of  the  nuclear  spins 
is  determined  by  minimizing  the  Helmholtz  free  energy, 

F  =  U  -TS,  (1.6) 

where  U  is  the  internal  energy  of  the  system,  T  the  temperature,  and  S  the  entropy. 

At  negative  temperatures,  the  situation  takes  a  peculiar  turn.  The  stability  crite- 
rion for  the  Helmohltz  free  energy  now  leads  to  the  system  maximizing  its  free  energy. 
This  can  be  seen  by  considering  the  thermodynamic  stability  of  the  system  shown  in 
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System  S 


Resevoir  S 


Composite 
System 


Figure  1.1:  The  composite  system  which  is  described  in  the  text. 

Figure  1.1.  The  composite  system  consists  of  a  large  thermal  reservoir,  and  the  sys- 
tem of  interest,  between  which  only  heat  can  be  exchanged.  The  reservoir  is  simply 
a  conceptual  tool  introduced  to  hold  the  temperature  constant.  The  total  entropy  of 
the  system  is  Stotai  =  Sr  +  S,  where  Sr  is  the  entropy  of  the  resevoir,  and  S  is  the 
entropy  of  the  system.  The  temperature  T  of  the  system  is  the  same  as  that  of  the 
resevoir,  which  is  at  Tr.  At  positive  temperatures,  the  equilibrium  of  the  composite 
system  is  characterized  by 


d{U  +  Ur)  =  0,  and 
d2(U  +  Ur)    >  0, 


(1.7) 
(1.8) 


given  that  dStotai  =  0,  which  implies  dSr  =  —  dS.  A  change  dSr  in  the  entropy  of  the 
resevoir  results  in  an  energy  change  dUr  =  Tr  dSr.  This  leads  to 


d(U  +  Ur)   =  0, 


=   dU  +  Tr  dSr, 
=  dU  +  T(-dS), 
=  d(U-TS), 


(1.9) 
(1.10) 
(1.11) 
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so  that  the  Helmholtz  free  energy  F  =  U  -  TS  of  the  system  is  an  extremal.  To 
determine  if  it  is  a  minimum  or  maximum,  first  consider  d2(U  +  Ur).  One  finds  that 

d2Ur     =     d^{dSr)\  (1-12) 

-  §-r(dSrf,  (1.13) 
=  0, 

as  the  reservoir  temperature  Tr  is  constant.  This  leads  to  d2(U  +  Ur)  =  d2U,  which 
is  greater  than  zero,  and  so  U  is  minimized.  The  condition  on  the  free  energy  now 
takes  the  form 

d2F   =   d2{U-TS),  (1.14) 
=  d2U-Td2S.  (1.15) 

The  differential  d2S  is  zero  since  S  is  the  independent  variable  whose  deviation  from 
equilibrium  is  being  considered.  This  leads  to 

d2F  =  d2U,  (1.16) 

so  that  d?F  is  greater  than  zero,  as  cPU  is  greater  than  zero.  This  implies  that  F  is 
a  local  minimum.  For  the  negative  temperature  case,  the  equilibrium  conditions  are 
set  by 

d(U  +  Ur)=0  (1.17) 

and 

d2(U  +  Ur)  <  0,  (1.18) 

at  constant  Stotai-  This  leads,  according  to  the  above,  to  the  free  energy  F  being  a 
local  maximum. 

However,  this  argument  concerning  the  free  energy  of  the  nuclear  spin  system  is 
not  quite  complete.  Based  on  the  above,  one  might  say  that  if  antiferromagnetic 


G 


(ferromagnetic)  order  is  achieved  at  positive  spin  temperatures,  then  ferromagnetic 
(antiferromagnetic)  order  is  the  natural  result  to  expect  for  ordering  at  negative 
temperatures.  This  holds  true  only  for  short  range  nearest-neighbor  interactions 
[13,  14].  The  long  range  dipolar  interaction  leads  to  the  possibility  of  having  both 
antiferromagnetic  or  ferromagnetic  ordering  at  opposite  signs  of  temperature.  This 
has  been  demonstrated  in  work  on  the  insulators  CaF2  and  LiF  by  Abragam  and 
co-workers  [13].  In  this  case,  the  negative  temperatures  were  achieved  in  the  rotating 
frame,  where  the  nuclear  spin  order  was  transferred  from  the  Zeeman  to  the  secular 
part  of  the  dipolar  Hamiltonian. 

Another  way  to  see  that  the  situation  is  not  straightforward  is  to  consider  the  inter- 
actions in  the  ordered  state  at  positive  and  negative  temperatures.  The  ordered  state 
at  negative  temperatures  is  identical  to  the  ordered  state  at  positive  temperatures  of 
the  system  in  which  the  signs  of  all  the  interactions  have  been  flipped.  Therefore,  if  a 
system  orders  ferromagnetically  at  positive  temperatures,  it  must  order  antiferromag- 
netically  at  negative  temperatures,  but  an  antiferromagnet  at  positive  temperatures 
can  be  an  antiferromagnet  at  negative  temperatures. 

A  phase  diagram  for  the  ordered  state  of  silver  nuclear  spins  at  negative  tem- 
peratures has  been  determined  experimentally  by  Hakonen  and  co-workers  [15].  The 
dashed  line  in  Figure  1.2  represents  a  theoretical  phase  boundary  for  negative  temper- 
atures, and  the  solid  line  is  the  phase  boundary  determined  from  [8]  for  the  antifer- 
romagnetically  ordered  state  at  positive  temperatures.  The  two  filled  circles  are  the 
basis  of  the  negative  temperature  phase  boundary  they  have  constructed.  However, 
the  paucity  of  data  makes  the  phase  diagram  speculative,  at  best.  From  their  data, 
they  estimate  the  critical  polarization  in  zero  field  and  5  //T  to  be  pcrit  —  —0.49, 
which  corresponds  to  a  critical  entropy  of  Scrit/R\n2  =  0.82.  The  two  data  points 
are  obtained  from  the  static  susceptibility  x'(0)  shown  in  Figure  1.3,  taken  in  zero 
magnetic  field  (Q),  and  at  5  ^T  parallel  (x)  and  perpendicular  (V)  to  the  sample 
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Figure  1.3:  Static  magnetic  susceptibility  in  zero  field  and  5  /xT  as  a  function  of  spin 
polarization,  from  reference  [14]. 
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foils.  The  static  susceptibility  is  calculated  from  the  absorption  signal  x"{f)  using 
the  Kramers-Kroonig  relation  x'(0)  =  (2/tt)  /0°°(x"(/)//)  df. 

The  polarization  scale  p  in  Figure  1.3  deserves  some  comment,  as  it  is  not  a 
meaningful  quantity  in  the  ordered  state,  or  in  zero  field.  The  polarization  of  the 
nuclear  spins  is  actually  determined  at  a  'high'  magnetic  field  of  200  fiT,  using  the 
formula  [13]  p  =  Af  x"(f)  df  (the  proportionality  constant  A  is  fixed  by  spectra 
taken  at  approximately  1  mK).  The  magnetic  field  is  then  adiabatically  reduced  to 
its  final  zero  (  or  low  )  field  value.  Instead  of  the  polarization,  the  horizontal  axis 
should  be  the  entropy  of  the  nuclear  spins,  which  is  well  defined  in  the  ordered  and 
paramagnetic  state.  For  non-interacting  spin  1/2  moments,  this  has  a  simple  relation 
to  the  absolute  value  of  the  polarization       given  by 

S/R  =  In 2  -  1/2((1  +  \p\)  ln(l  +  |p|)  +  (1  -  |p|)  ln(l  -  \p\)).  (1.19) 

Another  point  is  the  comparison  of  the  parameters  determined  from  the  Curie- 
Weiss  law 

where  C  is  the  Curie  constant,  and  0  is  the  Curie- Weiss  temperature.  Hakonen  and 
coworkers  found  0  =  —4.8  nK  at  positive  temperatures,  and  0  =  —2.8  nK  at  negative 
temperatures.  Within  the  mean  field  theory,  it  has  been  shown  [14]  that  one  finds  the 
same  0  for  both  signs  of  temperature.  This  result  is  in  agreement  with  the  nuclear 
susceptibility  measurements  on  rhodium  [16],  where  the  Curie- Weiss  temperature  was 
found  to  be  the  same  at  positive  and  negative  temperatures.  The  questions  regarding 
the  work  on  silver  by  Hakonen  and  co-workers  prompted  us  to  begin  studying  this 
system,  to  see  if  a  consistent  picture  could  be  developed. 
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1.3    Negative  Spin  Temperatures 

There  are  several  techniques  that  have  been  used  to  generate  negative  tempera- 
tures in  nuclear  spin  systems.  In  normal  NMR,  this  is  typically  done  with  a  180  degree 
pulse.  At  the  Larmor  frequency  jB0  (  where  B0  is  the  static  field  ),  the  excitation 
Bexc  is  applied  for  a  time  t  such  that  the  precession  angle  6  satisfies  9  =  ^Bexct  =  tt. 

Another  well  known  technique  is  that  of  adiabatic  fast  passage  [17].  A  magnetic 
field  B  of  fixed  magnitude  but  of  variable  direction  is  applied  to  the  magnetic  mo- 
ments. The  magnetization  M  is  initially  parallel  to  B.  The  changing  direction  of  the 
field  can  be  described  by  an  angular  velocity  u.  If  the  condition 

jB  >  u  (1.21) 

holds,  then  the  magnetization  M  will  stay  aligned  with  the  field  B,  even  as  it  moves. 

This  can  be  exploited  directly  in  high  frequency  pulsed  NMR  experiments.  A 
rotating  field  B\  of  frequency  u  is  applied  perpendicular  to  the  static  field  B.  When 
off-resonance,  the  magnetization  M  is  almost  completely  parallel  to  the  effective  field 

Beff  =  \JB\  +  [co/'j  —  Bo}2  in  the  rotating  frame.  By  sweeping  the  excitation  fre- 
quency through  the  resonance  slowly,  the  magnetization  M  will  follow  Beff  in  the 
rotating  frame.  At  resonance,  the  magnetization  is  aligned  with  B\.  By  continuing 
the  frequency  sweep  through  the  resonance,  the  magnetization  will  end  up  in  the  op- 
posite direction  of  the  static  field.  The  only  caveat  is  that  the  spin-lattice  relaxation 
time  (7\)  must  be  longer  than  the  time  it  takes  to  make  the  frequency  sweep. 

An  interesting  technique  was  used  by  Purcell  and  Pound  [18]  in  their  pioneering 
work.  They  removed  a  LiF  sample  from  a  large  static  field  and  placed  it  in  a  small 
solenoid  parallel  to  a  field  of  ~  100  gauss.  They  then  discharged  a  capacitor  through 
the  solenoid  so  that  the  direction  of  the  field  was  reversed  to  ~  -100  gauss,  and 
returned  the  sample  to  the  large  static  field  (  the  entire  process  took  approximately 
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3  seconds  ).  They  found  the  resulting  magnetization  to  be  reversed,  indicating  that 
the  moments  were  now  at  negative  temperatures. 

This  technique  was  refined  by  Oja  and  co-workers  [19],  and  has  subsequently  been 
used  as  the  chief  means  of  generating  negative  spin  temperatures  in  metals.  It  has 
been  employed  by  Hakonen  and  co-workers  in  their  work  on  silver  [12]  and  rhodium 
[16].  The  static  field  is  flipped  on  a  very  short  time  scale  compared  to  T2,  the  spin- 
spin  relaxation  time,  so  that  the  moments  do  not  follow  the  direction  of  field.  In 
Figure  1.4,  the  upper  part  of  the  diagram  shows  the  situation  before  the  field  is 
inverted,  and  the  lower  part  shows  the  final  result.  For  silver,  T2  ~  10  ms,  and  in 
rhodium,  T2  ~  10.5  ms.  The  best  inversion  efficiency  in  silver  was  90%  at  polarizations 
of  ~  30%,  while  at  higher  polarizations  (60-70%),  the  efficiency  dropped  to  ~  60-70%. 
In  rhodium,  for  polarizations  less  than  40%,  inversion  efficiencies  as  high  as  95%  could 
be  achieved.  However,  at  higher  polarizations,  the  inversion  efficiency  again  dropped 
below  60%. 

A  more  'sophisticated'  way  to  flip  the  spins  is  by  using  a  suitable  square  wave 
pulse.  This  is  essentially  a  single-period  version  of  the  180  degree  pulse,  and  was  first 
described  by  Friedmann,  et  al.  [20]  in  their  work  on  3He.  The  method  illustrates 
why  the  pulse  frequency  has  to  increase  as  the  pulse  duration  is  reduced  to  one  or 
two  periods  of  oscillation.  This  technique  was  also  tried  in  the  rhodium  work,  but 
efficiencies  of  only  60%  could  be  achieved,  even  at  low  polarizations  [16]. 

While  there  are  several  techniques  available  to  invert  the  nuclear  magnetization 
in  metals,  none  of  them  are  very  satisfactory  from  a  practical  standpoint.  Neither  the 
conventional  180  degree  pulse  nor  the  adiabatic  fast  passage  is  useful  for  bulk  metals 
at  very  low  temperatures.  They  require  high  frequencies,  which  leads  to  the  skin 
effect,  and  the  RF  field  Bexc  heats  up  the  conduction  electrons,  causing  the  nuclear 
spins  to  relax.  The  inversion  schemes  discussed  here  are  not  well  characterized,  and 
their  failings  at  high  polarization  are  not  understood.  The  inability  to  attain  high 
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Figure  1.4:  Schematic  representation  of  the  field  inversion  technique.  Note  that  the 
gyromagnetic  ratio  in  silver  is  negative.  Therefore,  the  spin  is  anti-parallel  to  the 
magnetic  moments,  which  are  depicted  with  arrows. 
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negative  polarizations  using  present  techniques  limits  the  work  that  can  be  done  on 
nuclear  spin  systems  in  the  ordered  state  at  negative  temperatures.  This  has  led  us  to 
develop  a  new  technique  for  inverting  the  nuclear  magnetization.  The  physical  model 
and  experimental  work  characterizing  this  new  scheme  are  discussed  in  Chapter  3. 
These  results  should  prove  useful  to  workers  in  the  low  temperature  field  of  nuclear 
magnetism  in  metals.  In  Chapter  2,  the  cryogenic  apparatus  and  electronics  used  in 
this  work  are  described.  Chapter  4  has  some  suggestions  for  future  work.  Appendix 
1  presents  the  results  of  measurements  on  the  nuclear  susceptibility  of  PrBe13,  and 
the  observation  of  nuclear  antiferromagnetic  ordering  in  this  system. 


CHAPTER  2 
EXPERIMENTAL  DESIGN  AND  TECHNIQUES 


To  cool  the  silver  nuclei  to  ultra-low  spin  temperatures,  the  silver  sample  was 
used  as  the  second  stage  of  a  two- stage  nuclear  demagnetization  cryostat.  In  the  first 
part  of  this  chapter,  the  cryogenic  apparatus  used  to  obtain  these  temperatures  is 
described.  The  second  section  covers  the  electronics  and  design  of  the  silver  experi- 
ment. The  third  section  describes  the  thermometry  used  to  monitor  the  performance 
of  the  cryostat.  The  final  section  discusses  the  experimental  procedures  used  in  the 
experiment. 

2.1  Cryostat 

The  cryostat  is  suspended  from  a  triangular  aluminum  table  that  rests  on  three 
Vibraplane1  air  isolation  pistons.  The  pistons  are  anchored  to  a  concrete  tripod  whose 
legs  are  isolated  from  the  building  foundation  and  rest  on  2  ton  concrete  pads.  The 
pads  serve  as  a  first  stage  in  de-coupling  the  cryostat  from  building  vibrations.  By 
pressurizing  the  pistons  with  N2  from  a  high  pressure  gas  cylinder,  one  effectively 
floats  the  table  and  mechanically  isolates  the  cryostat  from  the  tripod.  In  this  way, 
coupling  of  the  refrigerator  to  building  vibrations  is  further  minimized. 

Special  care  has  gone  into  providing  RF  shielding  of  the  cryostat  in  order  to 

minimize  noise-induced  heating.  A  large  screened  room  encloses  the  upper  part  of 

the  cryostat.  In  this  room  are  kept  all  of  the  sensitive  electronics  that  are  essential 

to  running  experiments  and  the  cryostat.  Computers  are  connected  to  the  equipment 
kinetic  Systems,  Boston,  MA  02131. 
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in  the  room  via  optical  cables  that  are  fed  through  a  small  pipe  in  the  roof  of  the 
shielded  room.  The  lower  part  of  the  cryostat  is  housed  in  a  large  aluminum  can  that 
is  bolted  to  the  floor  of  the  shielded  room. 

2.1.1    Dilution  Refrigerator 

The  cryostat  used  in  this  work  consists  of  two  principle  cooling  stages:  a  commer- 
cial dilution  refrigerator  unit  and  a  copper  nuclear  demagnetization  stage. 

The  dilution  refrigerator  is  an  Oxford  Instruments2  1000.  Dilution  refrigerators 
have  been  in  use  now  for  a  number  of  years,  and  the  general  principles  governing  their 
operation  are  well  documented  [21,  22],  so  only  those  aspects  key  to  this  work  will  be 
described.  A  more  thorough  description  of  the  cryostat  can  be  found  in  [23,  24,  25]. 

The  lowest  operating  temperature  of  the  dilution  unit  is  ~  4.5  mK  at  low  cir- 
culation rates  (~  0.35  mmol/s  of  3He).  This  base  temperature  is  typical  of  modern 
dilution  refrigerators.  However,  the  real  figure  of  merit  for  such  machines  is  their 
cooling  power.  In  Figure  2.1,  the  cooling  power  as  a  function  of  temperature  is  plot- 
ted for  a  fixed  3He  circulation  rate  of  h  =  1.27  mmol/s.  The  fit  is  according  to  the 
model  described  in  [26].  The  functional  form  of  the  cooling  power  can  be  written  as 

f  =  aT2-fc,  (2.1) 

where  a  and  b  are  fitting  parameters.  The  parameter  b  represents  the  residual  heat 

leak  to  the  fridge  at  this  circulation  rate,  and  is  found  to  be  1.7  fiW.  The  circulation 

rate  of  h  =  1.27  mmol/s  is  found  to  be  the  most  efficient  for  pre-cooling  the  nuclear 

demagnetization  stage.  This  upper  limit  for  the  pre-cooling  efficiency  occurs  because 

one  begins  to  circulate  appreciable  amounts  of  4 He  above  this  point,  while  the  3 He 

circulation  rate  becomes  saturated. 

2Oxford  Instruments,  130A  Baker  Avenue,  Concord,  MA  01742,  phone:  (508)  369-9933. 
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Figure  2.1:  Cooling  power  of  the  dilution  refrigerator. 

The  thermal  link  between  the  dilution  refrigerator  and  copper  nuclear  demag- 
netization stage  is  made  with  an  aluminum  superconducting  heat  switch  [27].  The 
superconducting  transition  temperature  for  aluminum  in  zero  magnetic  field  is 
TC(B  =  0)  ~  1.2  K,  while  the  critical  magnetic  field  at  0  K  is  BC(T  =  0)  ~  11  mT.  For 
temperatures  sufficiently  lower  than  the  transition  temperature  and  in  zero  applied 
field,  the  switch  is  thermally  insulating.  The  application  of  a  field  greater  than  Bc 
drives  the  aluminum  into  the  normal  state,  where  the  thermal  conductivity  due  to 
the  conduction  electrons  recovers  its  linear  temperature  dependence. 

2.1.2    Copper  Nuclear  Demagnetization  Stage 

The  demagnetization  stage  consists  of  nineteen  copper  sheets,  each  a  quarter-inch 
thick  and  36  inches  long.  The  width  of  each  plate  is  such  that  when  stacked  together, 
they  form  a  cylinder  2.5  inches  in  diameter.  The  top  part  of  the  plates  is  electron- 
beam  welded  to  a  large  gold-plated  copper  flange  4.5  inches  in  diameter,  which  is 
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bolted  into  a  copper  cage.  The  cage  is  designed  to  provide  space  for  experiments 
and  thermometers  to  be  attached  to  the  top  of  the  copper  bundle,  while  providing  a 
convenient  way  to  connect  the  bundle  to  the  bottom  of  the  aluminum  heat  switch. 

The  bottom  of  the  copper  bundle  ends  in  a  small  flange  to  which  the  silver  ex- 
periment is  bolted.  The  copper  flange  has  been  gold-plated  to  prevent  oxidation.  In 
Figure  2.2,  the  copper-to-silver  flange  connection  is  shown  in  detail.  Eight  3/4  inch 
6-32  Be-Cu  screws  are  used  to  hold  the  flanges  together.  To  ensure  good  metallic 
contact  at  low  temperatures,  small  tungsten  washers  (50  mils  thick)  were  used  to 
compensate  for  the  thermal  contraction  of  the  connecting  materials. 

The  end-compensated  8  T  magnets  that  are  used  for  adiabatic  demagnetization  of 
the  copper  bundle  and  silver  experiment  were  made  by  Cryomagnetics3.  The  magnet 
used  for  the  bundle  has  a  current  to  field  ratio  B/I  =  0.0935  T/A,  while  the  second 
magnet  for  the  silver  experiment  has  a  ratio  of  B/I  =  0.0976  T/A.  The  field  profile 
of  the  first  magnet  is  such  that  the  nuclear  heat  capacity  of  350  moles  of  copper  is 
equivalent  to  that  of  roughly  170  moles  in  the  peak  field  region. 

The  magnets  were  energized  to  75  amps  by  a  100  amp  Kepco  ATE6-100A  power 
supply4.  For  currents  less  than  10  amps,  a  smaller  Kepco  ATE  6-10M  unit  was  used. 
Both  power  supplies  were  operated  by  computer  controlled  DC  voltage  ramps  built  by 
the  physics  department  electronics  shop.  Details  on  the  demagnetization  procedure 
are  deferred  to  chapter  3,  where  experimental  procedures  are  discussed. 

2.2    Silver  Stage 

The  silver  experiment  is  in  essence  a  scaled  down  version  of  the  copper  stage.  It 
forms  a  second  demagnetization  stage  that  is  pre-cooled  by  the  copper  bundle  to  100 
//K  in  7  T.  In  Figure  2.2,  the  main  components  of  the  silver  stage  are  shown.  The 

3Cryomagnetics,  Inc.,  739  Emory  Valley  Road,  Oak  Ridge,  TN  37831-0548,  phone:  (615)  482- 
9551. 

4Kepco,  Inc.,  131-38  Sanford  Avenue,  Flushing,  NY  11352,  phone:  (718)  461-7000. 
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ure  2.2:  Silver  cross  and  sample  geometry. 
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Table  2.1:  Annealing  results  for  silver  foils. 


nT'\' 

KUJN 

T  (  C) 

Po2(Torr) 

time  (hours) 

nnn 

n  n  1\ 

Comments 

1 

650 

io-4 

5 

400 

turn  off  02  at  650  °C 

2 

700 

2 • IO"4 

4 

470 

turn  off  02  at  700  °C 

3 

750 

10~4 

3 

480 

turn  off  02  at  750  °C 

4 

800 

io~4 

1.7 

460 

turn  off  02  at  800  °C 

5 

800 

IO"4 

2.5 

675 

turn  off  02  at  room  temperature 

6 

800 

io-4 

2.5 

560 

turn  off  02  at  300  °C 

silver  flange  bolted  to  the  bottom  of  the  copper  stage  is  made  of  5  9's  purity  silver 
plate.  To  this  are  TIG- welded5  a  symmetric  pair  of  silver  stand-offs  that  clear  the 
copper  button  that  protrudes  from  the  bottom  flange  of  the  copper  stage.  The  silver 
cross  is  TIG-welded  to  these  standoffs.  The  cross  was  machined  from  a  5  9's  purity 
Marz  grade  silver  rod  from  MRC6.  The  cross  design  minimizes  the  amount  of  silver 
that  can  give  rise  to  eddy-current  heating  during  the  demagnetization  of  the  sample. 
It  also  provides  excellent  rigidity,  thereby  reducing  heating  due  to  vibrations.  The 
silver  flange  and  cross  were  annealed  for  10  hours  at  750  °C  under  an  oxygen  pressure 
of  3  •  10~4  torr.  From  previous  work  annealing  silver,  the  residual  resistivity  ratio 
(RRR)  of  the  components  is  estimated  to  be  >  550. 

The  bottom  of  the  silver  cross  ends  in  a  small  cup  into  which  the  sample  has  been 
welded  with  an  acetylene  torch.  The  silver  sample  consists  of  a  stack  of  16  silver  foils 
3  mm  x  50  mm  x  0.1  mm.  These  were  cut  from  a  100  mm  x  100  mm  x  0.1  mm 
Puratronic  grade  (  5  9's  pure  )  silver  foil  from  Alpha-Aesar7.  To  avoid  adding  any 
iron  impurities  to  the  sample,  a  pair  of  ceramic  scissors  from  Kyocera8  were  used  to 
cut  the  foils  to  the  proper  dimensions.  As  silver  has  a  high  vapor  pressure  at  elevated 

5  Tungsten  Inert  Gas 

"Materials  Research  Corporation,  560  Route  303  South,  Orangeburg,  NY  10962,  phone:  914-359- 
4200. 

7Aesar,  892  Lafayette  Rd.,  Seabrook,  NY  03874-1087. 

8Kyocera  America,  Inc.  8611  Balboz  Avenue,  San  Diego,  CA  92123. 
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temperatures,  some  tests  were  run  to  determine  the  optimum  annealing  procedure 
in  order  to  minimize  sample  loss  during  heat  treatment.  Table  2.1  summarizes  the 
results  of  this  work.  The  goal  was  to  maximize  the  residual  resistivity  ratio  while 
minimizing  the  annealing  time  to  cut  down  sample  loss.  The  sample  was  annealed 
using  the  parameters  from  run  5. 

To  minimize  the  effects  of  stray  magnetic  fields  on  the  experiment,  the  cold  plate 
radiation  shield  of  the  dilution  refrigerator  was  modified  to  seat  two  high  permeability 
cylinders  centered  on  the  silver  sample.  On  the  outside  is  fitted  one  of  mu-metal  ( 
0.1  in  wall  thickness,  7  in  long,  3.8  in  o.d.  )  and  on  the  inside  of  the  shield  is  one 
of  Cryoperm9  10  (  0.12  in  wall  thickness,  5  in  long,  3.63  in  o.d.  ).  To  keep  the 
permeability  of  the  material  high  at  low  temperatures,  the  Cryoperm  was  annealed 
in  >  1  atmosphere  of  4He  according  to  the  following  recipe  [22]:  1)  2  hours  at  1100 
°C,  2)  cool  to  500  °C  in  3  hours,  3)  cool  to  room  temperature  in  12  hours,  3)  anneal 
for  30  minutes  at  550  °C,  4)  cool  to  470  °C  in  40  minutes,  5)  anneal  for  2  hours  at 
470  °C,  6)  cool  to  room  temperature  in  2.5  hours. 

2.3    NMR  Coils  and  Electronics 

Continuous  Wave  (CW)  SQUID10  NMR  was  used  to  measure  the  magnetic  sus- 
ceptibility of  the  silver  nuclear  spins.  The  detection  system  consisted  of  two  sets  of 
orthogonal  coils  as  shown  in  Figure  2.3.  However,  the  first  RF  SQUID11  detector 
that  was  connected  to  the  vertical  pick-up  coil  never  worked  properly.  All  of  the  mea- 
surements reported  in  this  work  are  for  a  configuration  with  the  second  RF  SQUID 
connected  to  the  horizontal  pick-up  coil,  with  the  vertical  excitation  and  compensa- 
tion coils,  and  the  horizontal  static  field  coil.  Coil  parameters  are  listed  in  Table  2.2. 

9Vakuumschmelze,  D-6450,  Hanau,  Germany. 
10 Superconducting  QUantum  Interference  Device. 

nS.H.E  Corporation,  4174  Sorrento  Valley  Blvd.,  San  Diego,  CA  92121. 
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Figure  2.3:  The  relative  orientation  of  the  pick-up  and  excitation  coils. 
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Table  2.2:  Coil  parameters 


Design 

Wire 

B/I  (mT/A) 

Type 

Compensation  z 

4  mil  Cu 

1.3 

— 

Solenoid 

Compensation  x 

4  mil  Cu 

1.4 

— 

Helmholtz 

Excitation  z 

4  mil  Cu 

8.4 

Solenoid 

Excitation  y 

4  mil  Cu 

3.6 

Helmholtz 

Static  Field  z 

2  mil  Niomax 

40 

Solenoid 

Static  Field  y 

2  mil  Niomax 

3.55 

Helmholtz 

Pick-up  x 

2  mil  Niomax 

2.1 

Astatic  Helmholtz  pair 

Pick-up  z 

2  mil  Niomax 

1.98 

Astatic  solenoid  pair 

The  typical  experimental  arrangement  is  shown  in  Figure  2.4.  An  SRS12  DS340 
Digital  Signal  Generator  is  used  to  feed  both  the  excitation  and  compensation  coils. 
The  resistors  make  the  DS340  perform  as  a  current  source,  and  the  ratio  of  the  resistor 
values  is  chosen  to  balance  out  as  much  of  the  background  signal  as  possible.  The 
inverter  in  the  compensation  line  is  necessary  to  adjust  the  relative  phase  of  the  two 
signals  in  order  to  reduce  the  background.  The  pick-up  coils  form  a  superconducting 
loop  that  is  flux-linked  to  the  RF  SQUID  sensor.  The  SQUID  is  operated  in  the  slow 
flux-locked-loop  mode,  corresponding  to  a  maximum  slew  rate  of  1.3  x  103  (fosec, 
where  <j>0  is  the  flux  quantum.  The  feedback  signal  from  the  SQUID  electronics  is 
fed  into  an  Ithaco13  393  lock-in  that  uses  the  SYNC  signal  from  the  DS340  frequency 
generator  as  its  reference  input.  The  time  constant  setting  of  the  lock-in  is  1.25 
seconds  and  the  sensitivity  is  3  mV.  The  in-phase  and  quadrature  voltage  signals  from 
the  lock-in  were  measured  by  Hewlett-Packard  34401A  and  34574  digital  voltmeters 
(DVM).  A  PC  outside  of  the  shielded  room  using  LabWindows14  and  a  IEEE  488 

card15  is  connected  to  the  DVM's  and  signal  generator  via  an  HP-IB  optical  extender. 

12Stanford  Research  Systems,  Sunnyvale,  CA  94089. 

13Ithaco,  Inc.  735  W.  Clinton  St.,  Ithaca,  NY  14850. 

"National  Instruments,  6504  Bridge  Point  Parkway,  Austin,  TX  78730. 

15  ibid. 
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A  program  written  in  the  LabWindows  environment  controlled  the  frequency  and 
amplitude  of  the  signal  generator  while  recording  the  DVM  readings. 

2.3.1    Spin-flip  Electronics 

To  invert  the  silver  nuclear  moments,  a  large  sinusoidal  pulse  (comparable  to  the 
static  field)  is  sent  into  the  excitation  coil.  The  pulses  are  generated  with  the  Ar- 
bitrary Waveform  Composer  software  that  comes  with  the  Stanford  Research  DS345 
Digital  Signal  Generator.  A  PC  dedicated  to  running  this  software  package  is  con- 
nected to  the  DS345  via  another  HP-IB  optical  extender.  Once  the  desired  waveform 
has  been  generated,  it  is  loaded  into  the  DS345's  memory,  and  the  signal  genera- 
tor used  for  frequency  sweeps  is  disconnected  from  the  excitation  and  compensation 
coils.  The  'tipping'  generator  is  then  connected  to  the  excitation  coil  through  a  blQ 
resistance,  and  the  tipping  pulse  is  triggered  from  the  front  panel  of  the  DS345. 
The  tipping  generator  is  then  disconnected,  and  the  signal  generator  for  frequency 
sweeping  is  re-connected  to  the  compensation  and  excitation  coils. 

2.4  Thermometry 

The  temperature  of  the  copper  nuclear  demagnetization  stage  is  monitored  by  two 
thermometers.  A  3He  melting  curve  thermometer  (MCT)  is  used  for  temperatures 
above  1  mK,  while  a  platinum  (Pt)  NMR  thermometer  is  used  for  lower  temperatures. 

2.4.1    Melting  Curve  Thermometer 

The  melting  pressure  cell  is  a  Straty-Adams  strain  gauge  [28]  that  forms  one  of 
the  low  temperature  arms  of  a  capacitance  bridge.  The  electronic  configuration  of 
the  thermometer  is  shown  in  Figure  2.5. 
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Figure  2.5:  Electronic  circuit  for  the  melting  curve  thermometer. 


The  MCT  is  used  as  a  self-calibrating  primary  thermometer,  and  is  operated  in 
zero  magnetic  field.  The  fixed  points  used  to  calibrate  the  thermometer  are  the  super- 
fluid  A-phase  transition  (T^  =  2.505  mK)  [29]  in  liquid  3He  and  the  Neel  transition 
{TN  =  0.934  mK)  [29]  in  the  solid. 

The  strain  gauge  is  driven  by  the  internal  oscillator  of  an  Ithaco  393  lock-in  ampli- 
fier at  a  frequency  of  1.1  kHz,  with  an  amplitude  of  1.22  Vrms.  An  Eaton16  one-to-one 
isolation  transformer  is  used  to  separate  the  bridge  from  power  ground.  A  General 
Radio  1493  ratio  transformer  is  used  to  balance  the  bridge,  and  the  ratio  reading  is 
converted  to  temperature  based  on  a  capacitance  vs.  pressure  calibration  taken  at 
~  1  K.  The  small  off-balance  bridge  voltage  is  fed  to  the  lock-in,  and  the  resulting 
DC  output  is  sent  to  a  chart  recorder.  For  details  on  the  design  and  operation  of 
these  strain  gauges,  the  curious  reader  is  directed  to  the  review  article  by  Adams [30]. 


16Eaton  Corporation,  Electronic  Instrumentation  Division,  5340  Alia  Road,  Los  Angeles,  CA 
90066,  phone:  (213)  822-3061. 
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2.4.2    Pt  NMR  Thermometer 

The  Pt  thermometer  is  used  to  measure  the  temperature  of  the  copper  bundle 
below  the  Neel  transition  of  the  3He  solid,  where  the  slope  of  the  melting  curve 
tends  to  zero.  This  thermometer  has  no  fixed  points,  and  so  operates  as  a  secondary 
thermometer  calibrated  against  the  MCT  between  TV  and  ~  10  mK.  It  relies  on  the 
Curie  law  susceptibility  of  the  195Pt  spins,  and  can  be  used  to  measured  temperatures 
on  the  order  of  several  [iK. 

A  commercial  PLM-317  unit  was  used  to  perform  NMR  on  the  Pt  sample.  The 
layout  of  the  thermometer  and  electronics  is  shown  in  Figure  2.6.  The  cable  ca- 
pacitance of  the  coax  in  the  cryostat  plus  a  1  foot  section  of  BNC  cable  at  room 
temperature  gives  a  capacitance  to  ground  of  228  pF.  The  NMR  coil  that  completes 
the  tank  circuit  was  wound  with  3  mil  copper  wire  and  has  an  inductance  of  5.9  mH. 
The  tuning  capacitor  in  the  preamp  is  used  to  set  the  resonant  frequency  to  125  kHz. 
When  tuned,  the  tank  circuit  has  a  quality  factor  (Q)  ~  8.  The  static  field  coil  to  fix 
the  Larmor  frequency  of  the  Pt  nuclear  spins  is  made  from  9  mil  Formvar  insulated 
copper  clad  Supercon18  superconducting  wire.  The  magnet  consists  of  10  full  layers 
on  a  cyclindrical  copper  former  with  two  four  layer,  0.32  inch  long  end  correction 
windings.  A  minimal  amount  of  Stycast  1266  epoxy  was  used  to  hold  the  wires  in 
place.  The  magnet  sits  inside  a  copper  tube  whose  outer  surface  has  been  coated  with 
Pb,  which  is  superconducting  below  11  K.  The  entire  magnet  assembly  is  thermally 
anchored  to  the  mixing  chamber,  and  thermally  isolated  from  the  copper  bundle  by 
four  0.25  in  diameter  x  0.4  in  long  rock  maple  [31]  stand-offs.  The  current  to  fix  the 
Larmor  frequency  at  125  kHz  in  the  static  field  coil  is  1.26  A. 

It  was  found  that  the  measured  free-induction  decay  time  T2*  for  the  Pt  was 
only  0.6  ms,  which  compares  unfavorably  with  the  known  spin-spin  relaxation  time 

17RV-Elektroniikka  Oy,  Veromiehentie  14,  SF-01510  Vantaa,  Finland. 
18Supercon,  Inc.,  830  Boston  Turnpike,  Shrewsbury,  MA  01545. 
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Figure  2.6:  Schematic  diagram  of  the  Pt  NMR  circuit. 


of  T2  =  1.2  ms  [7].  This  was  attributed  to  field  inhomogeneities  in  the  magnet, 
which  were  reduced  by  inserting  a  superconducting  'jellyroll'  inside  the  static  field 
coil.  This  design  has  been  developed  by  Hechtfischer  [32],  and  exploits  the  boundary 
condition  of  the  magnetic  field  at  the  surface  of  a  flat  superconductor  to  produce 
very  homogeneous  small  fields.  The  jellyroll  is  made  of  a  3.5  inch  x  10  inch  x  2  mil 
sheet  of  NbTi  from  Teledyne  Wah  Chang19,  and  an  insulating  layer  of  1  mil  thick 
Mylar.  The  Mylar  and  NbTi  are  wound  around  a  copper  cylinder  several  times,  and 
care  is  taken  so  that  there  are  no  superconducting  shorts  between  the  NbTi  layers. 
Small  amounts  of  Stycast  1266  and  dental  floss  were  used  to  hold  the  jellyroll  on  the 
former.  After  this  was  inserted  into  the  bore  of  the  static  field  coil,  T2*  was  found  to 
be  ~  1  ms. 

The  500  Pt  wires  that  form  the  NMR  thermometer  are  in  the  form  of  a  brush 
that  is  butt-welded  to  a  small  5  9's  pure  Pt  plate.  The  Pt  assembly  was  annealed 
in  air  for  3  hours  at  900  °C.  In  order  to  avoid  quartz  contaminating  the  sample,  the 


19 


Teledyne  Wah  Chang,  Albany,  OR  97321. 
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Figure  2.7:  Temperature  of  the  bundle  during  the  precool. 

Pt  was  placed  in  a  small  alumina20  boat  inside  a  quartz  tube  in  the  annealing  oven. 
This  recipe  yielded  RRR  values  of  over  1000  on  Pt  test  samples  [33]. 

2.5    Experimental  Procedures 


A  typical  run  to  take  NMR  data  on  the  silver  experiment  starts  with  precooling 
the  copper  bundle  in  a  field  of  approximately  7  T  for  3-4  days.  This  brings  the  tem- 
perature of  the  bundle  to  just  over  11  mK,  and  the  magnetic  heat  capacity  of  the 
copper  nuclei  is  on  the  order  of  200  J/K.  The  temperature  of  the  bundle  during  this 
procedure  is  shown  in  Figure  2.7.  The  bundle  is  then  demagnetized  from  these  initial 
conditions  according  to  the  following  recipe  distilled  from  practical  experience  and 
computer  simulations: 

y  rp  210mjn  g  rj  r-p  180mm  j  g  "J1  120mm  q  gg  rj-, 


20 


SGP,  Willow  Grove,  PA  19090. 
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When  the  field  of  the  magnet  is  ~  3  T,  the  heat  switch  is  opened  (made  superconduct- 
ing) and  the  circulation  rate  of  the  dilution  refrigerator  is  reduced  to  ~  0.35  mmol/s. 
The  magnet  is  persisted  at  0.93  T  in  order  to  switch  over  from  the  100  A  to  the  10  A 
Kepco  power  supply,  which  is  then  used  to  control  the  remainder  of  the  demagneti- 
zation for  the  copper  bundle  with  a  higher  precision. 

The  100  A  Kepco  is  now  brought  inside  the  RF  shielded  room  and  connected 
to  the  second  magnet.  The  demag  of  the  first  magnet  is  continued  with  the  10  A 
supply  to  0.37  T  in  150  minutes,  at  which  point  the  magnet  is  again  put  in  persistent 
mode  and  a  4He  transfer  is  made.  The  second  magnet  is  then  energized  to  6.7  T  and 
persisted.  The  final  stage  of  the  bundle  demagnetization  is  continued  according  to: 
0  37  T  120m^"  0  19  T  240m^n  q  094  T  240m^"  o  047  T 

The  copper  bundle  is  now  at  95  (J.K.  There  is  now  a  12  hour  wait  as  the  silver 
nuclei  come  into  thermal  equilibrium  with  the  copper  stage  and  reach  nearly  100% 
polarization. 

For  small  polarization  p,  the  relaxation  of  the  silver  nuclei  is  described  by  the  well 
known  rate  equation 

where  Ti  is  the  nuclear  spin-lattice  relaxation  time  and  peq  corresponds  to  the  equi- 
librium polarization.  Using  the  Korringa  relation  k  =  T\Teq  [34],  and  the  fact  that 
for  small  polarization  p  oc  1/T,  the  above  equation  can  be  re-written  as 

—  (  — )  =  Zk(  J_  _  —\  (2.3) 

dt    Teq  K     Teq  T 

where  Teq  denotes  the  temperature  to  which  the  silver  nuclear  spins  are  relaxing.  This 
is  identified  with  the  temperature  of  the  silver  conduction  electrons,  which  are  taken 
to  be  in  equilibrium  with  the  conduction  electrons  of  the  copper  bundle. 
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For  high  polarization,  the  above  rate  equation  no  longer  holds,  and  has  to  be 
modified  [35,  36]  to  take  into  account  the  full  relation  of  the  polarization  to  the 
magnetic  field  B.  The  rate  equation  then  becomes 

dp     h^B/2kB  (  ,,x 


dt  k 


(coth(|)-coth(^))p,  (2.4) 


where  x  =  h^B/hsT,  xeq  =  h^B /kBTeq,  and  /  is  the  nuclear  spin.  The  above 
equation  can  be  simplified  as  the  polarization  for  spin  I  =  1/2  is  given  by 

Substituting  this  back  into  the  rate  equation  yields 

dt  k       \p  peqJ 

In  Figure  2.8,  the  polarization  as  a  function  of  time  is  shown  for  both  the  107 Ag 
and  109 Ag  nuclear  spins  (time  zero  corresponds  to  the  demag  of  the  copper  bundle 
from  0.37  T).  Twelve  hours  after  the  demagnetization  of  the  copper  bundle,  the 
1  K  pot  is  closed  and  the  magnet  on  the  silver  experiment  is  demagnetized  from 
6.7  T  l2^n  0.04  T.  At  this  point,  a  small  2-A/2-V  power  supply  (  powered  by  deep 
cycle  marine  batteries  )  is  used  to  energize  the  small  NMR  static  field  coil  to  ~  90  //T, 
and  the  current  in  the  large  second  magnet  is  ramped  to  zero  in  1  minute.  The  static 
field  coil  is  then  increased  to  ~  850  //T.  The  large  magnet  is  persisted,  and  the  current 
leads  from  the  100  A  Kepco  are  disconnected.  The  large  magnet  is  now  depersisted 
with  the  current  leads  open  circuited  for  1  minute,  and  then  persisted.  At  this  point, 
the  flux-transformer  heater  for  the  SQUID  is  turned  off,  and  hopefully  everything 
is  working.  The  data  acquisition  program  is  now  initiated,  and  frequency  sweeping 
begins. 
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Figure  2.8:  Predicted  pre-cooling  curve  for  the  silver  nuclei  at  100  fjK  in  6.7  T. 


For  the  frequency  sweeps,  the  excitation  field  Bexc  is  typically  5  nT,  and  the 
canonical  frequency  step  is  2  Hz.  The  program  has  a  user  defined  dwell  time  for 
each  step,  which  is  set  to  2  or  3  seconds.  Once  the  dwell  time  has  been  reached,  the 
program  queries  the  DVM's  to  obtain  both  channels  of  the  lock-in  output,  and  saves 
the  information  to  a  file.  The  program  then  steps  the  frequency  to  the  next  value. 
After  acquiring  the  raw  data,  there  is  usually  an  appreciable  amount  of  signal  left. 
To  eliminate  this,  a  number  of  large  tipping  pulses  are  used  to  heat  the  spins  until 
the  magnetization  is  in  the  noise.  At  this  point,  frequency  sweeps  of  the  background 
signal  are  taken.  Typically,  four  background  sweeps  are  acquired  for  signal  averaging. 
The  excitation  signal  is  now  disconnected,  and  several  sweeps  with  the  compensation 
signal  alone  are  taken.  This  gives  a  'calibration'  of  the  frequency  dependence  of  the 
gain  of  the  SQUID,  as  the  compensation  coil  does  not  couple  to  the  sample,  but  does 
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couple  to  the  SQUID  pick-up  coil.  Note  that  all  of  the  data  acquired  is  taken  in  the 
same  static  field,  as  the  background  has  a  large  field  dependence. 

After  the  data  taking  is  complete,  the  1  K  pot  is  opened  (  the  maximum  'hold' 
time  of  the  pot  is  ~  7  hours  ).  The  flux-transformer  heater  is  turned  on,  and  the 
small  static  field  coil  is  ramped  to  zero.  The  copper  bundle  is  then  remagnetized 
according  to: 


A  helium  transfer  is  made  at  the  completion  of  the  remagnetization,  and  the  second 
magnet  is  energized  to  6.7  T  and  persisted.  The  demag  of  the  copper  bundle  is  then 
carried  out  from  0.37  T  as  before,  and  the  whole  process  cycles  through.  A  single 
demagnetization  of  the  large  copper  bundle  is  capable  of  precooling  the  silver  stage 
4-5  times  in  this  fashion. 


CHAPTER  3 

THEORETICAL  CONSIDERATIONS  AND  EXPERIMENTAL  RESULTS 


The  steady-state  response  of  the  silver  nuclei  are  modelled  using  the  Bloch  equa- 
tions. The  coupling  of  the  two  silver  nuclear  species  via  the  RK  interaction  and 
the  dipolar  interaction  for  the  sample  foil  geometry  are  included.  Lineshapes  for  the 
sample  are  calculated  for  the  particular  excitation  and  pick-up  coil  assemblies  used 
in  this  work.  The  theoretical  description  of  the  spin-flip  dynamics  is  then  described, 
and  the  last  sections  compare  the  experimental  results  with  the  models. 

3.1    Modified  Bloch  Equation 

The  magnetic  properties  of  nuclear  moments  are  well  described  by  the  phenomono- 
logical  equations  of  Bloch  [37].  In  the  case  of  highly  polarized  silver,  where  there  are 
two  distinct  magnetic  nuclear  species,  the  equations  have  to  incorporate  an  extra  term 
that  describes  the  RK  and  dipolar  coupling  between  the  two.  This  term  is  defined  as 

R  =  E^V  (3-D 

PoP7  n    i  rij 

and  is  the  dimensionless  measure  of  the  strength  of  the  exchange  interaction  relative  to 
the  averaged  dipolar  interaction.  Here  is  the  strength  of  the  exchange  interaction, 
and  p  is  the  number  density  of  silver  atoms.  The  equations  of  motion  for  the  107 Ag 
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and     Ag  magnetizations  now  take  the  form  (ignoring  relaxation  terms  for  simplicity) 


dM 


107 


dt 


=  7107  MW7  x 


B0  +  Bexc  +  (R+  1/3)  /i0M109+  D  Mo(M107  +  M109) 


(3.3) 


and 


dM 


109 


dt 


7109^09  X 


B0  +  Bexc  +  (R  +  1/3)  tioM107+  D  /i0(M107  +  Mm) 


,  (3.4) 


where  D  is  the  shape-dependent  demagnetization  factor  [17],  the  1/3  arises  due  to 

— *  — * 

the  Lorentzian  field,  and  B0  and  Bexc  are  the  static  and  excitation  fields,  respectively. 
This  form  of  the  Bloch  equations  was  first  used  by  Ekstrom  and  coworkers  in  their 
work  on  copper  [1]. 

The  outline  of  the  calculations  will  be  for  a  sample  with  cylindrical  geometry,  as 
this  is  analytically  easier  to  handle,  and  the  physical  meaning  of  the  results  is  clearer. 
The  modifications  of  the  solutions  to  the  experimental  foil  geometry  will  be  described 
at  the  end. 


3.1.1    Stcadv-State  Solutions  for  the  Wire  Geometry 


For  the  wire  geometry,  the  demagnetization  factor  is  diagonal,  and  has  Dxx  = 
Dyy  —  — 1/2  and  Dzz  =  0.  The  linear  excitation  field  Bexc  is  decomposed  into  counter- 
rotating  fields,  and  the  only  significant  response  will  be  from  the  counter-clockwise 
component  Bexc  —  Bi(cos(u>t),  —  sin(cj^),  0). 

The  condition  for  steady-state  precessions  is  given  by  approximating  the  i-component 
of  the  107 Ag  and  109 Ag  magnetizations  with  their  respective  equilibrium  values,  so 
that  M2107  ^  M0m  and  M]m  w  M0109.   For  ease  of  calculation,  complex,  instead 
of  vector  notation,  is  now  used.   The  excitation  field  now  takes  the  form  B^xc  = 
Biexp(jivt).  The  precessing  magnetizations  are  given  by  M|07  =  m107  exp(jujt)  and 

M109  =  m109exp(^)}  where  M107,109  =  MK)7,109  +  jjtflOT.lW 
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The  equations  of  motion  now  take  the  form 


109 


jcomm   =   -jml07^07  [Bo  +  (R+  l/3)/i0M0109  +  (l/2)//0M0107 
+jmmll07(R  -  l/6VoM0107  +;M01077l07JB1 

jumim   =    -jm1097l09  [Bq  +  {R+  1/3)/x0M0107  +  (l/2)/x0M, 
+jm107ll09(R  -  l/6)/i0M0109  +  jMo109^^!. 

From  this  system  of  equations,  one  finds  eigenfrequencies 


(3.5) 


(3.6) 


u±  = 


2  L 


^io7  +  ^109  ±  ^Ko9  -  wi07)2  +  47io77io9(i?  -  1/6)2^M0107M0109 


,  (3-7) 


where 


and 


w107  —  7107 


B0  +  (i?  +  l/3)/z0M0109  +  (l/2)/x0M01 


rl07 


u/109  =  7l09  [b0  +  (R  +  l/3)/i0M0107  +  (l/2)/ioM0109] 


(3.8) 


(3.9) 


To  first  order,  the  resonance  frequencies  are  modified  in  two  ways.  In  silver, 
R  ~  —2.5  [3],  so  that  the  term  (R  +  l/3)//oM)  corresponds  to  an  overall  downshift 
due  to  the  antiferromagnetic  exchange  between  the  unlike  nuclear  spin  species.  The 
last  term,  (l/2)//0M)  gives  an  upshift  to  the  spins  due  to  the  torque  of  the  precessing 
demagnetization  field  of  the  same  spin  species.  This  is  illustrated  by  the  data  in 
Figure  3.1. 

To  higher  orders,  there  is  an  overall  repulsion  term  47l077io9(-R—  l/6)2//oM0107M0109 
that  represents  the  coupling  between  the  precessing  107 Ag  and  109 Ag  magnetizations. 
This  term  scales  as  (R  -  1/6)2M0107M0109.  When  the  polarization  p  of  the  107 Ag  and 
109 Ag  spins  are  the  same,  this  term  becomes  proportional  to  (R-  l/6)2p2,  the  square 
of  the  polarization  p,  as  M0  oc  p.  So  for  small  polarizations,  it  is  essentially  neglible. 
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Figure  3.1:  A  frequency  sweep  at  B  =  825.4  /zT  showing  the  downshift  and  suppres- 
sion/enhancement of  the  silver  signal  as  described  in  the  text.  The  solid  lines  mark 
the  unshifted  resonance  frequencies. 


The  solutions  for  the  steady-state  magnetizations  with  no  relaxation  take  the  form 

(3.10) 


m  107  /  r  „  mW  uioo  ~  w  +  7109 {R  ~  l/6)//0A/0109 
™    fo\  =  7107  M0   


((J  —  U!+){u>  —  UI-) 


and 


H 109  ^107  ~  ^  +  7iQ7(^  -  1/6)^0M0107 
(cj  —  w+)(u;  —  w_) 


m107£i  =  7iw  Af* 


(3.11) 


The  intensities  of  the  two  modes  shown  in  Figure  3.1  are  not  equal,  nor  do  they  scale 
as  the  ratio  of  the  nuclear  moments  of  the  107 Ag  and  109 Ag  spins.  This  result  can  be 
understood  by  studying  the  response  functions  found  above  at  resonance. 
At  the  lower  mode       for  the  107 Ag  spins  one  has 


m107/B1  =  -lmMt 


107  ^109  ~  fcJ-  +  7ios(fl  ~  l/gj/xo-Mp1 

}  (u+  -  W_) 


109 


-8(u}  —  UJ-), 


(3.12) 


where  S(x)  is  the  delta  function. 
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The  numerator  can  be  expanded  as 


=  1/2 


wiw  -  ^07+VVi09  -  oj[07Y  +  47io77io9(fl  -  1/6)2^M0107M0109|(3.13) 


+7io9(i?-l/6)/ioM109, 

~    ^109  ~  W107  ^  ?  ?  

^109  w107 

+7lo9(i?-l/6)/xoM109,  (3.15) 
>  0. 

This  last  expression  is  greater  than  zero,  so  that  the  overall  response  of  the  107 Ag 
nuclear  spins  at  u;_  is  negative.  Performing  a  similar  calculation  for  the  109 Ag  spins, 
one  finds  that  their  response  is  positive.  Extending  this  to  the  ui+  mode,  the  107 Ag  and 
109 Ag  responses  are  found  to  be  negative.  The  result  is  outlined  in  Figure  3.2.  The 
important  point  is  that  the  higher  frequency  mode  is  now  identified  with  an  in-phase 
precession  of  both  the  107  Ag  and  109  Ag  nuclear  spins,  while  the  lower  frequency  mode 
is  shown  to  be  where  the  109  Ag  nuclear  spins  precess  out  of  phase  with  respect  to 
the  107 Ag  spins.  The  sum  of  the  two  responses  reproduces  the  features  of  Figure  3.1, 
where  there  is  an  overall  enhancement  of  the  signal  intensity  of  the  upper  mode,  and 
a  suppression  of  the  intensity  of  the  lower  mode. 

3.1.2    Foil  Geometry 


The  experiment  in  this  work  has  a  sample  with  a  foil  geometry  and  uses  a  linear 
excitation.  The  demagnetization  factor  has  only  one  non-zero  component,  Dxx  =  —  1, 
which  breaks  rotational  symmetry.  The  linear  excitation  field  is  now  given  by  Bexc  = 
Bi(0,0,sin(ut)). 

In  the  case  of  the  wire  geometry,  the  resulting  magnetizations  are  circularly  po- 
larized. This  arises  due  to  the  neglect  of  the  counter-rotating  component  of  the 
excitation  field  in  converting  the  linear  excitation  to  a  rotating  one.  The  complete 
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Figure  3.2:  Relative  signs  of  the  107 Ag  and  109 Ag  nuclear  spins  at  the  resonance 
frequencies. 

solution  would  have  an  elliptically  precessing  magnetization  even  for  the  rotating 
excitation.  The  inclusion  of  relaxation  terms  further  complicates  matters,  as  two  dif- 
ferent phase  shifts  are  then  required  to  complelely  describe  the  magnetizations  in  the 
x  and  z  directions.  With  these  considerations  in  mind,  the  steady-solutions  for  the 
foil  geometry  take  the  form 

M™7  =  mmcos(ujt  +  a),  (3.16) 

Ml07  =  rW7m107sm(ujt  +  (3),  (3.17) 

Mz109  =  m109cos(^  +  7),  (3.18) 

M]m  =  r109mW9sm(ut  +  5),  (3.19) 

where  the  r's  and  phase  shifts  a,  {3,  7,  and  S  describe  the  elliptical  response  of  the 
magnetizations.   Substituting  these  solutions  into  the  Bloch  equations  for  the  foil 
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geometry,  one  obtains  the  substantial  system  of  equations 


1  /Ti  07 

0 

() 

0 

x 

0  " 

77^107  COs(of) 

1  / 7107 

(J 

0 

0 

0 

0 

13  r 

mio7  sin(o:) 

-Az 

0 

a; 

1/T107 

Bz 

0 

0 

0 

^107^107  C0S(/?) 

0 

-1/T107 

0 

Bz 

0 

0 

rio7?7iio7  sin(/?) 

0 

0 

Cx 

0 

a; 

1/^109 

"At 

0 

miog  cos(7) 

t) 

0 

0 

-1/7"109 

u 

0 

-Dx 

mio9  sin(7) 

cz 

0 

0 

0 

-A 

0 

1/T109 

no9"aio9  cos(^) 

0 

cz 

0 

0 

0 

-A 

-1/^109 

a; 

.  rio9m109  sin(5) 

7io7M0107 
0 
0 
0 

7io9M0109 
0 
0 
0 


The  various  terms  are 

Ax 
Az 
Bx 
Bz 
Cx 

cz 

Dx 


7l07  (B0  +  (^+l/3)/i0M0109), 

7107  (B0  +  (R+  l/3)/i0M0109  +  //0M0107), 

7107  (i?+l/3)/i0M0107, 

7107  (R  -  2/3)fi0MQ107, 

7io9(fl+l/3)//0M0109, 

7109  (i?-2/3)/i0M0109, 

7109  {B0  +  (R  +  l/3)/i0M0107),  and 

7109  (B0  +  (R  +  l/3)/i0M0107  +  //0M109). 


(3.20) 


As  the  pick-up  coil  in  this  experiment  is  parallel  to  the  i-axis,  the  detected  signal 
is  given  by 


Mf07  4-  Mf09   =   mio7  cos(cut  +  a)  +  m109  cos(ut  +  7) 
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=   —[^107  sin(of)  +  m109  sin(7)]  sin(u;t) 

+  [mi07  cos(q;)  4-  miog  cos(7)]  cos(o;t).  (3-21) 

While  it  is  not  practical  to  write  out  the  full  closed  form  solution  to  the  above 
system  of  equations,  for  the  case  of  no  relaxation  (r  — >  oo),  a  similar  analysis  as  for 
the  wire  geometry  can  be  carried  out.  The  ratio  of  the  signal  intensities  for  the  109 Ag 
and  107 Ag  modes  is  found  to  be 

7io7M0107Q1+  +  7i09M0109ft2+ 

r  7i07M0107J71_  +  7io9M010^2_  {6-ZZ) 

where 

=  oj2.  —  DXDZ  —  AZCX  —  DXCZ  —  BZCX, 

tt2+   =  —  AXAZ  -  AXBZ  -  DZBX  -  BXCZ, 

Qi_   =  u}_  —  DXDZ  —  AZCX  —  DXCZ  —  BZCX, 

f^2-   =  w_  —  AXAZ  —  AXBZ  —  DZBX  —  BXCZ, 


and  the  eigenfrequencies  for  the  foil  geometry  take  the  form 

«4   =   \[(AXAZ  +  DXDZ)  +  2BXCZ±((AXAZ-DXDZ)2  +  4[(AXAZ  +  DXDZ)BXCZ 

+AXDXBZCZ  +  AZDZBXCX})1/2].  (3.23) 

The  ratio  of  the  line  intensities  is  a  function  of  the  static  field  B0  and  the  p.  Given 
B0  and  the  ratio,  the  polarization  can  be  determined. 

3.2    Spin-flip  Equations  of  Motion 

In  conventional  pulsed  NMR,  the  magnetization  can  be  inverted  with  a  180°  pulse 
with  an  excitation  frequency  uexc  equal  to  the  Larmor  frequency  u0.  However,  for 
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extremely  short  pulses  (one  or  two  periods  long)  where  the  excitation  field  Bexc  is 
on  the  order  of  the  static  field  B0,  this  is  no  longer  the  case.  The  experiment  by 
Friedmann  and  co-workers  [20]  illustrates  that  for  single-period  pulses,  the  tipping 
frequency  should  be  greater  than  the  Larmor  frequency.  For  square- wave  pulses,  this 
is  uexc  =  \/2cj0,  as  they  have  shown. 

In  silver,  the  situation  is  complicated  owing  to  the  presence  and  interaction  of 
the  two  isotopes  107 Ag  and  109 Ag  .  At  low  fields,  where  the  two  resonance  lines  have 
exchange  merged  [3]  into  one,  even  a  single  period  pulse  turns  out  to  be  too  long 
in  comparison  to  T2  to  efficiently  flip  the  magnetization.  The  calculation  described 
in  this  section  maps  out  the  inversion  efficiency  as  a  function  of  tipping  field  and 
frequency  in  silver.  The  optimal  efficiency  is  found  to  occur  for  tipping  frequencies 
greater  than  the  average  Larmor  frequency  of  the  two  lines,  in  qualitative  agreement 
with  the  results  of  Friedmann.  However,  it  is  very  close  to  the  frequency  of  the  109 Ag 
mode,  due  to  the  coupling  between  the  two  isotopes. 

The  full  equations  of  motion  for  a  finite  tipping  pulse  are  described  by  the  Bloch 
equations.  A  suitable  change  of  coordinate  system  allows  the  inclusion  of  relaxation 
terms  in  the  description.  These  are  important  as  T2  ~  10  ms  in  silver,  and  the 
experimental  frequency  range  of  this  work  is  only  1  to  1.5  kHz,  at  which  even  a 
single-period  tipping  pulse  lasts  for  ~  1  ms. 

The  starting  point  for  the  description  of  the  tipping  spin  dynamics  is  the  Bloch 
equations.  However,  unlike  the  usual  case  of  pulsed  NMR,  where  a  transformation 
[6]  is  made  to  the  rotating  frame,  we  keep  the  full  set  of  equations  and  solve  them 
numerically.  The  usual  trick  describes  the  tipping  field  as  a  rotating  magnetic  field 
when  it  is  actually  a  linear  field.  This  transformation  ignores  the  counter-rotating 
component  of  the  tipping  field,  which  is  important  in  this  case  [38]. 
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3.2.1    Relaxation  Terms 

The  magnetization  relaxes  due  to  spin-lattice  (Ti)  and  spin-spin  (T2)  interactions. 
For  the  experiments  in  this  work,  7\  is  effectively  infinite  due  to  the  very  low  experi- 
mental temperatures,  and  can  be  ignored.  T2,  on  the  other  hand,  is  approximately  10 
ms.  As  the  experimental  frequencies  are  on  the  order  of  /  ~  1  kHz,  and  1/T2  ~  16  Hz, 
the  relaxation  of  the  magnetization  from  T2  processes  has  to  be  included  in  the  equa- 
tions of  motion  to  correctly  describe  the  spins  during  a  tipping  pulse.  It  is  easiest 
to  consider  only  one  spin  species  to  determine  how  the  spin-spin  relaxation  enters  in 
the  finite  tipping  pulse  case.  Then,  the  method  outlined  will  be  applied  to  the  case 
in  silver,  yielding  the  full  coupled  equations  of  motion. 

The  total  applied  magnetic  field  is  the  sum  of  the  DC  and  excitation  fields 

Btotal  =  B0  +  Bexc,  (3-24) 

where  B0  =  (0,0,  B0)  and  Bexc  =  (B\  sin(o>i),  0, 0).  In  the  standard  transformation 
[13,  17]  from  the  lab  frame  to  the  rotating  frame,  B0  is  replaced  by 

Bo+cuf-yz.  (3.25) 

For  most  NMR  work,  the  excitation  field  is  considered  to  be  a  superposition  of  two 
counter-rotating  fields  of  magnitude  B^/2  with  frequencies  ±o».  At  resonance,  when 
lo  ~  7.B0,  the  component  at  frequency  —  u  will  be  off  resonance  by  2u,  and  can  be 
ignored.  However,  when  Bx  ~  B0,  this  approximation  breaks  down.  To  account  for 
this,  a  transformation  of  the  Bloch  equation  must  be  made  to  the  'rocking'  frame  ( 
shown  in  Figure  3.3  ).  This  introduces  another  effective  field  term  .  It  also  allows  the 
incorporation  of  the  T2  relaxation  in  a  straightforward  manner.  The  equations  in  the 
rocking  (primed)  frame  now  take  the  form 


Mz,   =  Mx,6 


(3.26) 
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Figure  3.3:  Relation  between  the  effective  and  static  field.  The  Y  and  Y'-axes  are  the 
same,  and  point  into  the  page. 


Mx,   =   -Mz,6  +  jMy,Btotal 


M„ 


,   ,        „  My' 

-~/Mx,Btotal  -  -f- 


(3.27) 
(3.28) 


From  the  field  geometry  in  Figure  3.3,  one  has  that 


tan  8  = 


Bi  sm(ojt) 

Bn 


(3.29) 


Taking  the  time  derivative  of  both  sides  yields 


6  = 


Bi/Bqlu  cos(cj^) 


l  +  (Bl/B0)2sm2(ut) 
The  total  field  Btotal  is  now  given  by 


(3.30) 


Btotal  =  yjBi  +  Bf  sin2 (tot). 

The  transformation  back  to  the  lab  frame  is  effected  by 

Mz    -    Mz>  cos  9  -  Mx>  sin  0, 
sin  6  +  Mxi  cos  6, 

My    =  Myl, 


(3.31) 

(3.32) 
(3.33) 
(3.34) 
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where  the  initial  conditions  are 

Mz>  =  1,  (3.35) 

Mx,  =  0,  (3.36) 

My,  =  0.  (3.37) 

The  coupled  equations  of  motion  for  the  silver  nuclei  can  now  be  derived,  using  the 
previous  equations  as  a  guide.  After  some  algebra,  one  finds  that 

M2109   =   (cos^sin^/T2  +  7io9^o(i?+l/3)My107)Mr 
-7l09(Bi  sm(ut)  +tto(R+  l/3)j\4107)My109 
-  sin2  e/T2M]m.  (3.38) 

Mx109   =   -  cos2  6/T2M1xQ9+  7109CB0  +  M#  +  l/3)M2107)My109 

+(cos0  sin0/T2  -  7l09/io(^+  l/3)My107)Mz109,  (3.39) 

<9  =  -7109(^0+ MR+i/mr'wr-Mi09/^ 

+7l09(B1  sinful)  +  Ho(R  +  l/3)Mi07)M,109,  (3.40) 

As  the  original  set  of  equations  is  symmetric  with  respect  to  the  two  nuclear  spin 
species,  the  results  for  107 Ag  can  be  obtained  by  making  the  substitution  107 Ag 
<->  109  Ag  . 

These  equations  do  not  allow  a  closed  form  solution,  and  can  only  be  numeri- 
cally integrated.  A  program  using  The  Student  Edition  of  MATLAB  IV1  is  used 
to  perform  the  integration  for  different  values  of  the  parameters.  In  the  simulation, 
the  z-component  of  the  107 Ag  and  109 Ag  magnetizations  are  used  to  characterize  the 
efficiency  of  the  tipping  pulse.  The  total  Ml°tal  =  M2107  +  M2109  is  compared  before 
and  after  the  pulse,  giving  a  calculated  efficiency  of 

Mtotai  after  tipping  pulse  . 
calc     Mlotal  before  tipping  pulse ' 


iMathWorks,  Inc.,  24  Prime  Park  Way,  Natick,  MA  01760-1500,  Phone:  (508)  647-7000 
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The  results  of  the  simulation  show  several  key  points.  At  low  fields,  where  the 
Larmor  frequency  is  ~  100  Hz,  the  tipping  efficiency  for  one  and  two  period  pulses 
is  no  better  than  50%.  This  limitation  occurs  owing  to  the  tipping  frequency  being 
comparable  to  T2,  the  spin-spin  relaxation  time.  The  pulse  duration  is  long  enough  to 
allow  considerable  de-phasing  to  occur  among  the  spins.  The  only  way  to  increase  the 
efficiency  according  to  the  simulation  is  to  move  to  higher  fields  (  and  correspondingly 
higher  tipping  frequencies  ).  At  Larmor  frequencies  of  ~  1.5  kHz,  the  tipping  efficiency 
from  the  simulation  is  found  to  be  much  higher,  on  the  order  of  80%.  However,  the 
optimum  tipping  frequency  is  not  simply  related  to  the  resonance  frequencies  of  the 
107 Ag  and  109 Ag  spins  for  polarizations  of  at  least  10%.  To  achieve  a  high  efficiency 
inversion  of  the  magnetizations  of  both  silver  isotopes,  the  tipping  frequency  must 
almost  be  as  large  as  the  resonance  frequency  of  the  109  Ag  mode. 

3.3    Data  Reduction 

The  output  from  the  lock-in  has  been  de-convoluted  according  to  the  following 
scheme: 

1.  The  lock-in  offset  is  subtracted  from  the  raw  data,  background  data,  and  cali- 
bration data. 

2.  A  gain  correction  factor  of  2.887  is  made  to  the  cosine  channel  of  the  lock-in 
output,  to  compensate  for  extra  amplification  in  the  sine  channel. 

3.  The  phase  4>  and  amplitude  A  of  the  raw  data,  background  data,  and  compen- 
sation data  are  calculated  from  <j)  =  arctan(y/x)  and  A  =  y/x2  4-  y2  at  each 
frequency,  where  x  and  y  are  the  cosine  and  sine  channels  of  the  lock-in. 
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Figure  3.4:  The  lock-in  output  before  data  reduction.  The  cosine  and  sine  channels 
of  the  lock-in  are  denoted  by  +  and  Q>  respectively. 


4.  The  phase  and  amplitude  of  the  raw  and  background  data  are  unmixed  with 
the  calibration  data.  This  involves  calculating  new  phase  angles 

<f>'taw     =    draw  —  (Peat  ~  ^offset, 
^background     =     ^background  —  ^cal  ~  ^offseti 

where  <p0ffset  is  a  phase  correction  of  0.51  radians.  Both  phases  of  the  raw  data 
are  now  reconstructed  from 


x  =  {Araw/Acat)  x  cos  $ 


and 


(3.42) 


(3.43) 


y  =  (Araw/Acat)  X  Sin  <f>'rav). 

The  background  signal  is  reconstructed  in  the  same  manner. 
5.  The  background  is  then  subtracted  from  the  raw  data. 

A  set  of  raw  data  is  shown  in  Figure  3.4.  The  result  of  following  the  above  steps 
is  shown  in  Figure  3.5. 


/ 
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Figure  3.5:  Processed  data,  as  described  in  the  text. 

3.3.1    Field  Calibration 

The  static  field  coil  in  the  y  is  calibrated  by  scaling  the  107 Ag  and  109 Ag  peak 
frequencies  at  low  polarization  (  ~  0.05  )  against  their  published  gyromagnetic  ratios. 
From  [7],  one  finds 

7iQ7  1-732 
7109  1-991 

=  0.8699,  (3.44) 

where  the  ratio  of  the  Knight-shifted  values  has  been  used.  The  peak  frequencies  are 
determined  by  fitting  a  double  Lorentzian  to  both  phases  of  the  signal.  The  spectrum 
used  for  the  field  calibration  is  shown  in  Figure  3.6.  The  solid  line  is  the  result  of  a 
non-linear  least  squares  fit  to  both  phases  using  the  Nelder-Mead  simplex  algorithm 
[39].  The  ratio  of  the  fitted  resonance  frequencies  is 

1443 

/o109         1  659 
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Figure  3.6:  Data  used  to  determine  the  horizontal  static  field  calibration.  Solid  lines 
are  the  result  of  the  fit  described  in  the  text. 


=  0.8698 


(3.45) 


which  compares  well  with  the  accepted  value.  This  gives  a  field-to-current  ratio  of 
B/I  =  3.55  mT/A.  The  calculated  value  from  Table  2.2  is  3.55  mT/A. 

The  spin-spin  relaxation  times  from  the  Lorentzian  fit  are  5.5  ms  and  4.9  ms  for 
107 Ag  and  109 Ag  ,  respectively.  These  are  smaller  than  the  published  [40]  values  of  9.1 
ms  and  10.4  ms.  The  discrepancy  is  accounted  for  by  a  field  inhomogeneity  of  ~  1%, 
determined  from 

1  1 


To*  T, 


(3.46) 


where  T2*  is  the  experimentally  determined  spin-spin  relaxation  time,  T2  is  the  true 
spin-spin  relaxation  time,  and  AB  is  the  spread  in  the  static  field. 


48 


0.4 


0.2 


0.1 


a 
3 


>> 

13  _o.i 


a 

y 
3 

CO 


-0.2 


-0.3 


-0.4 


 1  1  1  1  1  1  1  

    f  \  

\  - --"V- 

IS                 +  N- 

r  *** 

i  i  i  i  1  

'1200        1250    .    1300        1350        1400        1450        1500       1550  1600 

Frequency  (  Hz  ) 

Figure  3.7:  Spectrum  used  to  set  the  polarization  calibration.  The  solid  line  is  the  fit 
described  in  the  text. 


3.3.2    Polarization  Calibration 


The  processed  spectra  are  all  fit  to  double  Lorentzians  as  described  in  the  previous 
section.  The  fitted  lineshapes  for  the  107 Ag  and  109 Ag  modes  are  then  integrated,  and 
the  intensity  ratio  of  the  two  modes  is  used  to  determine  the  polarization  as  described 
at  the  end  of  section  3.1.2. 

An  attempt  was  made  to  measure  the  equilibrium  polarization  at  100  fjK.  This 
would  have  provided  an  'absolute'  measure  of  the  polarization  to  the  integrated  sig- 
nal intensity.  Unfortunately,  the  noise  level  is  presently  too  large  to  determine  this 
number.  The  polarization  scale  in  this  work  has  been  set  by  selecting  one  spectrum 
and  performing  the  above  analysis  to  fix  the  polarization.  The  spectrum  used  to  fix 
the  polarization  scale  is  shown  in  Figure  3.7.  The  solid  lines  are  the  fit  to  the  data, 
which  have  been  integrated  to  give  the  ratio  of  the  line  intensities.  The  measured 
lineshape  falls  off  faster  than  the  fit,  indicating  that  the  functional  form  is  closer 
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Figure  3.8:  Initial  magnetizations  at  positive  temperature  for  the  tipping  sequence. 
Solid  line  is  a  Lorentzian  fit  to  the  data,  which  gives  p  =  0.27. 

to  a  Gaussian,  rather  than  Lorentzian,  lineshape.  The  integration  was  carried  out 
numerically  using  the  trapezoid  technique  [41].  The  static  field  for  this  data  set  is 
777.4  pT,  and  the  signal  intensity  analysis  yields  a  polarization  p  =  0.27.  This  gives 
a  calibration  constant  of 


V 


=  10.2  x  10 


-3 


area 


(3.47) 


3.4    Spin-flip  Results 


In  Figure  3.8,  a  spectrum  taken  at  B0  —  832.5  pT  is  shown.  Figure  3.9  shows  the 
magnetization  after  the  application  of  a  tipping  pulse  of  frequency  1625  Hz  applied 
for  two  periods  with  a  vertical  excitation  field  of  Bexc  =  540  pT  (  from  the  calculated 
coil  constant  in  Table  2.2  ).  The  initial  polarization  is  pi  =  0.27.  and  is  taken  to  be 
the  same  for  both  spin  species.  After  tipping,  the  polarization  from  the  integrated 
area  is  found  to  be  pj  =  —0.22,  giving  an  efficiency  of  ~  81%.  The  magnitude  of  the 
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Figure  3.10:  Final  magnetizations  at  positive  temperature  for  the  tipping  sequence. 
Solid  line  is  a  Lorentzian  fit  to  the  data,  which  gives  p  =  0.16. 
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tipping  field  is  used  as  a  fitting  parameter  to  match  the  initial  and  final  polarizations 
in  the  simulation,  as  there  is  no  calibration  of  the  excitation  coil.  Using  a  tipping  field 
of  Bexc  =  450  //T  in  the  simulation,  the  final  calculated  polarization  for  the  107 Ag 
and  109 Ag  spins  is  found  to  be  pcatc  =  —0.21,  yielding  a  model  efficiency  of  78%. 

The  magnetization  is  now  flipped  from  negative  temperature  back  to  positve  tem- 
perature by  applying  the  same  pulse  as  before.  The  resulting  spectrum  is  shown  in 
Figure  3.10.  The  final  calculated  polarization  for  the  107 Ag  and  109 Ag  spins  is  found  to 
be  slightly  different,  giving  pw7  =  0.16  and  p10g  =  0.14.  The  polarization  determined 
from  the  integrated  intensity  is  found  to  be  p  =  0.16. 


CHAPTER  4 
FUTURE  WORK  AND  CONCLUSIONS 

4.1    Higher  Initial  Polarization 

The  initial  polarizations  for  the  nuclear  magnetization  reported  in  this  work  are 
only  of  moderate  value,  being  in  the  area  of  ~  0.2-0.3.  From  the  discussion  in  Chapter 
I,  this  is  too  low  a  starting  point  to  reach  the  ordered  state  of  the  nuclear  spins  at 
negative  temperatures,  for  which  the  zero-field  critical  polarization  [15]  is  reported  as 
pcrit  =  0.49. 

There  are  several  modifications  that  can  be  tried  to  determine  the  cause  of  these 
low  initial  polarizations.  In  the  current  detection  scheme,  the  final  Larmor  field  ap- 
plied to  the  sample  is  horizontal,  while  the  field  from  the  large  second  demagnetization 
magnet  is  vertical.  This  leads  to  the  possibility  of  polarization  loss  due  to  inadequate 
control  in  the  changing  of  the  magnetic  field  orientation  from  vertical  to  horizontal 
at  the  end  of  the  demag.  To  try  and  go  smoothly  from  one  orientation  to  the  other, 
the  small  Larmor  field  coil  in  the  vertical  direction  should  be  energized  along  with 
the  horizontal  field  coil,  towards  the  end  of  the  second  stage  demag.  The  second 
stage  magnet  should  then  be  ramped  down  in  the  normal  way,  after  which  the  small 
vertical  field  can  be  reduced  smoothly. 

Another  possibility  is  to  swap  the  vertical  excitation  and  horizontal  static  field 
coils.  That  is,  use  the  horizontal  static  coil  for  the  excitation  field,  and  the  vertical 
excitation  coil  for  the  static  field,  while  keeping  the  same  pick-up  coil  orientation.  This 
is  something  that  could  be  done  rather  easily  from  the  top  of  the  cryostat  at  room 
temperature.  The  only  difficulties  might  arise  from  the  large  diameter-to-length  ratio 
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of  the  vertical  coil  (  about  5  to  1  ),  giving  poor  field  homogeneity  and  the  possiblility 
of  a  large  phase  shift  in  the  excitation  from  the  horizontal  static  coil.  One  ancillary 
benefit  from  the  swapping  of  coils  is  that  it  could  enable  the  calibration  of  the  current- 
to-field  ratio  of  the  excitation  coil.  This  would  remove  one  fitting  parameter  from  the 
tipping  simulation. 

4.2  Conclusions 

In  this  work,  a  new  technique  has  been  developed  to  generate  high  efficiency  inver- 
sions of  the  nuclear  magnetization  in  metals.  The  experimental  results  in  silver  have 
been  found  to  be  in  good  agreement  with  the  model  developed  from  the  Bloch  equa- 
tions using  a  linear  excitation  field  and  incorporating  the  RK  exchange  interaction 
between  the  nuclear  moments.  The  best  experimental  inversion  efficiency  achieved  in 
this  work  is  ~  81%  for  a  starting  polarization  p  =  0.27,  yielding  a  final  polarization  of 
p  =  —0.22.  This  technique  has  not  yet  been  fully  optimized,  as  the  simulation  shows 
that  the  inversion  efficiency  increases  at  higher  static  fields  for  a  suitable  choice  of 
tipping  frequency  and  magnitude.  With  this  new  scheme,  it  should  be  possible  to 
reach  the  ordered  state  of  the  silver  nuclear  moments  at  negative  temperatures,  and 
determine  the  phase  diagram. 


APPENDIX  A 
SUSCEPTIBILITY  MEASUREMENT  OF  PRBE 


This  appendix  is  a  brief  summary  of  work  done  on  the  hyperfine-enhanced  nuclear 
paramagnet  PrBei3.  The  first  part  of  this  section  describes  the  experimental  appa- 
ratus used  in  the  susceptibility  measurements,  along  with  the  results.  The  next  part 
gives  a  brief  introduction  to  the  physics  of  this  material,  and  presents  the  results  of 
a  mean  field  calculation  used  to  interpret  the  data. 

A.l  Introduction 

Hyperfine  enhanced  Van  Vleck  paramagnets  have  been  a  subject  of  interest  to 
the  low  temperature  community  for  a  number  of  years.  Their  use  as  possible  nuclear 
refrigerants  [42]  for  low  temperature  cryostats,  along  with  the  interesting  physics  [43] 
in  these  systems,  has  made  them  an  ongoing  field  of  study. 

Recent  work  by  Hiess  and  co-workers  [47]  on  the  NpBe13  system  has  highlighted 
the  possibility  of  understanding  the  heavy-fermion  (HF)  superconductor  UBe13  by 
studying  the  other  rare-earth  and  actinide  Metal-Bei3  compounds.  These  systems 
have  a  cubic  NaZnJ3  structure,  where  the  Metal  component  occupies  the  sites  of  a 
simple  cubic  lattice.  All  of  the  Uranium  based  HF  compounds  are  found  to  order 
antiferromagnetically,  with  TN  >  Tc,  with  the  glaring  exception  being  UBe13  [48,  49]. 
That  is,  except  for  UBei3,  the  other  Uranium  based  HF's  all  exhibit  strong  magnetic 
correlations.  The  neutron  scattering  work  of  Hiess  and  co-workers  finds  a  magnetic 
wave  vector  q  =  (1/3  0  0)  for  NpBei3,  the  same  as  found  in  several  other  of  the 
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Figure  A.l:  Silver  platform  and  coils  for  the  PrBe13  experiment. 

Metal-Be13  systems  [47].  They  suggest  that  this  is  the  wave  vector  that  should  be 
searched  for  magnetic  correlations  in  UBe13. 

This  leads  to  the  interesting  possibility  of  understanding  the  magnetism  in  UBe13 
through  studying  the  magnetic  structure  of  PrBei3.  Of  all  the  Metal-Bei3  compounds, 
PrBei3  is  unique,  as  it  is  the  4/  analog  of  UBei3.  The  electronic  configuration  of  U4+ 
is  5/2  (  assuming  it  is  not  in  a  mixed- valence  state  ),  while  Pr3+  has  a  configuration  of 
4/2.  Neutron  scattering  experiments  on  PrBe13  should  provide  interesting  information 
relevant  to  the  physics  of  UBe13. 

Previous  work  on  PrBe13  has  been  focused  on  high  temperatures  [44,  45],  where 
only  the  ionic  moments  of  the  system  have  been  studied,  or  to  methods  [46]  where 
the  interaction  between  the  Pr  nuclear  moments  was  indirectly  determined.  The 
experiment  described  here  directly  probes  the  Pr  nuclear  spin  interaction. 
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A.  1.1    Sample  and  Experimental  Apparatus 

The  zero-field  nuclear  magnetic  susceptibility  of  PrBei3  was  measured  with  a 
mutual  inductance  technique  [50]  using  a  DC  SQUID  sensor.  The  polycrystalline 
sample  weighs  3.17  mg,  and  has  a  volume  of  ~  1  mm3.  The  sample  was  provided  by 
Professor  G.  R.  Stewart  of  the  University  of  Florida. 

A  stand  (  see  Figure  A.l  )  consisting  of  four  bunched  silver  wires  welded  to  a 
silver  plate  formed  the  thermal  link  and  sample  holder  for  the  PrBe13.  The  sample 
was  glued  to  the  wires  using  Epo-Tek1  conductive  silver  paint.  The  stand  was  bolted 
to  the  top  of  the  copper  nuclear  demagnetization  cryostat,  with  the  wires  and  sample 
inserted  into  the  excitation  and  detection  coils. 

The  experimental  coil  arrangement  is  outlined  in  Figure  A.l.  The  sample  is  cen- 
tered in  the  top  of  an  astatic  pair  of  pick-up  coils.  As  the  wire  forms  a  closed 
superconducting  loop,  flux  is  conserved  on  the  detection  side  of  the  circuit.  The  total 
flux  in  this  part  can  be  written  as 

<f>=  (MPl(l  +  X)  +  MP2)ip  +  MC2ic,  (A.l) 

where  Mpi  is  the  mutual  inductance  between  the  primary  coil  and  the  first  astatic 
coil,  MP2  is  the  mutual  inductance  between  the  primary  coil  and  second  astatic  coil, 
MC2  is  the  mutual  inductance  between  the  compensation  coil  and  second  astatic  coil, 
and  ip  and  ic  are  the  currents  in  the  primary  and  secondary  coils.  The  mutual 
inductance  between  the  primary  and  first  astatic  pair  is  modified  by  the  presence 
of  the  magnetic  sample  by  a  factor  of  (1  +  x(T)),  where  x(T)  is  the  temperature 
dependent  susceptibility  of  the  nuclear  moments.  Denoting  the  inductances  in  the 
detection  part  of  the  circuit  as  LUL2  and  L3,  the  total  flux  can  also  be  written  as 

4>  =  (Lx  +  L2  +  L3)is,  (A.2) 
JEpoxy  Technology  Inc.,  14  Fortune  Dr.,  Billerica,  MA  01821. 
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Figure  A. 2:  Schematic  of  the  electronics  for  the  PrBe13  experiment. 
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where  ia  is  the  super-current  flowing  in  the  circuit.  These  equations  can  be  combined 
to  find 

x       Sis  T,  U 

5x  =  1 — iu~'  A-3 
ip  MPl 

8is  is  related  to  the  DC  SQUID  output  voltage,  which  the  manufacturer  gives  as 
1.832  xlO-5  A/V.  The  mutual  inductance  MpX  is  3.3  //H,  and  the  total  inductance  £  Li 
is  4.9  IjlH.  The  primary  current  is  just  ip  =  Vexc/Rp,  where  Vexc  is  the  excitation  voltage 
from  the  HP  synthesizer,  and  Rp  is  the  current  limiting  resistor  to  the  excitation  coil. 
The  excitation  voltage  was  set  at  0.6  Vpp.  The  excitation  field  used  in  this  work  has 
an  amplitude  of  2.64  ^Tpp  at  a  frequency  of  9.23  Hz.  With  this  in  hand,  one  can 
convert  the  lock-in  voltage  to  a  susceptibility. 

Figure  A. 3  shows  the  results  of  the  measurement.  The  susceptiblity  has  a  cusp  at 
TV  =  0.40  mK,  indicating  that  the  Pr  nuclear  moments  have  ordered  antiferrmagnet- 
ically.  The  high  temperature  part  of  the  nuclear  susceptibility  data  was  fitted  to  a 
Curie- Weiss  form,  plus  a  constant  background  term.  A  non-linear  least  squares  fit  to 
the  data  at  high  temperature  yields  a  Weiss  temperature  0  =  —  2.2  mK.  The  results 
of  the  fit  are  shown  in  Figure  A. 3.  This  measured  Weiss  temperature  is  in  contra- 
diction to  the  ESR  g-shift  result  [46],  which  found  a  weak  ferromagnetic  interaction 
corresponding  to  0  =  60  //K. 

A.2    Mean  Field  Calculation  for  Tw 


In  the  mean  field  theory,  the  ratio  of  the  Weiss  to  the  Neel  temperature,  Q/TN, 
can  be  used  to  determine  the  arrangement  of  the  moments  in  the  ordered  state.  The 
structure  that  our  experimental  result  implies  has  been  calculated  using  a  mean  field 
theory  incorporating  nearest-neighbor  (nn)  and  next-nearest-neighbor  (nnn)  interac- 
tions on  a  simple  cubic  lattice.  The  procedure  used  follows  that  outlined  by  Smart 
[51],  and  is  in  agreement  with  the  result  of  ter  Haar  [52]. 
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Figure  A. 3:  Nuclear  susceptibility  of  PrBei3.  Solid  line  is  the  non-linear  least  squares 
fit  to  the  high  temperature  data  from  5  to  20  mK.  The  ordering  transition  is  shown 
in  the  inset. 
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Figure  A. 4:  Division  of  the  cubic  lattice  into  eight  sub-lattices,  as  described  in  the 
text. 

The  Heisenberg  Hamiltonian  for  the  nuclear  spins  is  of  the  form 

H  =  -^JijIiIj,  (A.4) 

i<j 

where  is  the  exchange  interaction  between  spins  and  the  S  are  spin  angular 
momentum  operators.  The  lattice  is  divided  up  into  n  =  8  sub-lattices  relative  to  the 
center  spin  1,  such  that  no  interactions  take  place  between  any  neighbors  on  their  own 
sub-lattice.  The  division  of  spins  is  shown  in  Figure  A.4,  where  the  numbers  denote 
the  sites  belonging  to  each  sub-lattice.  Note  that  the  sub-lattice  corresponding  to 
number  8  is  for  next-next-nearest-neighbors,  for  which  the  interactions  are  zero.  It  is 
introduced  only  for  convenience  [51]  in  the  following  calculations. 
The  total  magnetic  field  acting  on  the  $h  sub-lattice  is  given  by 

n 

Hi  =  HQ  +  Y,'YiA,  (A.5) 

where  H0  is  the  external  field,  Mj  is  the  magnetization  of  the  fh  sub-lattice,  and  7„ 
is  the  coupling  between  the  sub-lattices.  All  self-coupling  terms  of  the  form  ju  are 
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taken  to  be  zero.  In  terms  of  the  Hamiltonian,  one  has 

7»j  =  n(zijJij/{NA{gnnn)2),  (A.6) 

where  Zy  is  the  number  of  j  neighbors  of  spin  i,  jj,n  is  the  nuclear  magneton,  Na  is 
Avogadro's  number,  gn  is  the  Lande  g-factor,  and  n  is  the  number  of  sub-lattices.  The 
usual  mean  field  coefficients  are  related  to  these  7's  by  7^  =  /jfc7i,  where  /jfc  is  the 
fraction  of  the  #h  nn's  of  spin  j  on  sub-lattice  A;.  The  normal  mean  field  coefficients 
are  found  to  be 

7nn  =  3znnJnn/(NA(gnHn)2),  (A.7) 

where  znn  is  the  number  of  nearest  neighbors  and  Jnn  is  the  interaction  strength 
between  nearest  neighbors.  The  form  of  jnnn  is  obtained  by  replacing  nn  with  nnn 
in  the  above  equation. 

By  putting  the  external  field  Hq  equal  to  zero  in  equation  A. 5,  one  obtains  a 
set  of  equations  relating  the  sub-lattice  magnetizations  and  their  local  fields.  The 
eigenvalues  of  this  set  of  equations  yield  the  possible  mean  field  ordering  temperatures 
for  the  magnetic  moments. 

The  ordering  temperatures  as  a  function  of  the  mean  field  coefficients  are  found 
to  be 

T\  =  {jrinn  ~  lnn)C ,  (A.8) 
T2  =  -l/S{jnn  +  Jnnn)C,  (A.9) 

T3  =  1/3(7™  -  1nnn)C,       and  (A.10) 

T4  =  (7nn  +  7rmn)C,  (A.ll) 

where  C  =  NAHo(9nVn)2I(I  +  l)/{3kB)  is  the  nuclear  Curie  constant.  The  Weiss 
temperature  is  0  =  (7„n  +  7nnn)C,  corresponding  to  the  temperature  T4.  Ordering 
temperatures  one,  two,  and  three  correspond  to  different  possible  antiferromagnetic 
structures.  The  structure  corresponding  to  T3  can  be  eliminated  as  it  describes  an 
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Figure  A. 5:  The  ratio  of  the  Weiss  and  Neel  temperatures  as  a  function  of  the  mean 
field  coefficients  and  the  spin  configurations  for. 
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energetically  unfavorable  arrangement  of  spins,  as  compared  to  the  Neel  temperatures 
T\  and  Ti- 

In  Figure  A. 5,  the  ratio  of  the  Weiss  temperature  0  to  the  Neel  temperature  TN 
is  plotted  as  a  function  of  the  mean  field  coefficients  7,  and  the  dashed  and  solid  lines 
correspond  to  7\  and  T2,  respectively.  The  lowest  energy  configurations  for  the  spins 
correspond  to  the  upper  branches  of  the  plot.  The  results  of  the  calculation  show  that 
the  ground  state  for  the  nuclear  spins  is  of  the  Type  II  structure  for  ynnn/7nn  >  1/2, 
while  for  jnnn/lnn  <  1/2  the  ordered  structure  is  Type  I.  From  the  data,  one  finds 
that  Q/TN  ~  —6,  which  implies  antiferromagnetic  interactions  between  nnn  and  a 
Type  II  ordered  structure  for  the  Pr3+  nuclear  moments.  The  large  value  (  by  about 
a  factor  of  2  )  of  the  experimentally  determined  ratio  as  compared  to  the  result  of 
the  calculation  suggests  that  the  transition  is  not  mean  field  like,  and  fluctuations  are 
important. 


APPENDIX  B 
CW  NMR  MEASUREMENTS  ON  GOLD 


Two  different  high  purity  gold  samples  were  studied.  For  both  samples,  no  NMR 
signal  was  detected.  The  following  sections  describe  the  experimental  geometry  of 
the  gold  samples,  and  the  heat  treatments  used  to  prepare  them.  The  electronics  and 
NMR  coils  used  for  the  gold  experiments  are  the  same  as  those  described  in  chapter 
2. 

B.l    Gold  Fan  Stage 

The  thermal  link  between  the  bottom  of  the  copper  demagnetization  stage  and 
this  sample  consisted  of  a  four  9's  purity  gold  sheet1  that  was  folded  over  seven  times 
into  the  shape  of  a  fan.  The  width  of  the  folds  was  0.25  inches,  and  the  overall 
length  of  the  fan  was  6.86  inches.  The  five  9's  purity  gold  sample2  consisted  of  16 
foils  0.1  mm  x  3  mm  x  50  mm.  These  were  cut  from  a  50  mm  x  50  mm  x  0.1  mm 
sheet  with  a  pair  of  ceramic  scissors.  The  foils  were  diffusion  welded  into  the  bottom 
end  of  the  fan,  with  the  overlap  of  the  sample  and  fan  being  10  mm.  The  top  part 
of  the  gold  fan  was  diffusion  welded  to  a  set  of  nine  silver  fins  (  1.125  in  x  0.2  in  x 
0.016  in  )  that  were  TIG- welded  into  the  silver  base  attached  to  the  bottom  of  the 
copper  demagnetization  stage.  The  overlap  between  the  silver  fins  and  the  gold  fan 
was  1  in. 

The  gold  fan,  gold  sample,  and  silver  base  were  all  annealed  separately  before  dif- 
fusion welding.  The  silver  stage  was  annealed  for  10  hours  at  a  temperature  of  750  °C 

'A.D.  MACKAY,  Inc.,  10  North  Broadway,  P.O.  Box  G,  Red  Hook  NY  12571. 
2Aesar,  892  Lafayette  Road,  RO.  Box  1087,  Seabrook,  NH  03874-1990. 
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Table  B.l:  Annealing  results  for  gold  foils. 


RUN 

T  (°C) 

Po2(Torr) 

time 

RRR 

Comments 

1 

900 

26  hours 

100 

annealed  in  air 

2 

925 

IO"4 

21  hours 

140 

3 

965 

io-4 

7  days 

365 

4a 

965 

10"4 

10.25  days 

445 

4b 

965 

10"4 

16.6  days 

488 

4c 

965 

io-4 

20.1  days 

445 

cumulative  time  on  sample  4 

under  an  oxygen  pressure  of  3  x  10-4  torr.  The  results  of  various  heat  treatments  on 
gold  foils  are  summarized  in  Table  B.l.  The  gold  sample  foils  were  annealed  for  9 
days  at  950  °C  under  an  oxygen  pressure  of  10-4  torr.  The  gold  fan  was  annealed  in 
the  same  fashion,  except  that  small  quartz  slides  were  placed  in  the  folds  of  the  fan 
to  keep  it  from  welding  itself  together  during  the  heat  treatment. 

The  sample  foils  were  then  diffusion  welded  to  the  end  of  the  gold  fan.  A  small 
stainless  steel  clamp  with  quartz  slide  'jaws'  was  used  to  hold  the  sample  foils  and 
fan  in  place.  This  heat  treatment  was  done  at  900  °C  for  ten  hours  under  vacuum 
(~  4  x  10~6  torr  ).  At  some  point  during  the  diffusion  welding,  a  leak  opened  up  in 
the  system,  resulting  in  an  oxide  layer  (  from  the  stainless  steel  clamp  )  forming  on 
the  gold  sample.  This  was  removed  using  a  mild  solution  of  HBr.  The  gold  fan  was 
then  diffusion  welded  to  the  silver  fins  by  cooking  at  660  °C  for  12  hours  under  an 
oxygen  pressure  of  10-4  torr.  After  cleaning,  the  gold  foils  were  clean  enough  to  stick 
to  each  other  on  contact.  In  order  to  keep  the  foils  separate,  a  silicon  oxide  powder  ( 
six  9's  pure  )  from  Alpa-Aesar  was  sprinkled  on  the  foils.  The  powder  adhered  lightly 
to  the  gold,  but  was  sufficient  to  keep  the  foils  from  bonding  to  each  other. 
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B.2    NIST  gold  sample 

The  second  gold  experiment  used  the  same  design  of  silver  thermal  link  as  in  the 
silver  experiment.  The  geometry  of  both  experiments  is  identical,  the  only  difference 
being  that  the  silver  sample  foils  in  Figure  2.2  are  replaced  by  gold  foils. 

A  small,  high  purity  (  six  9's  )  gold  rod  was  purchased  from  NIST3  (  Standard 
Reference  Material  685  ).  At  first,  several  attempts  were  made  to  roll  the  gold  into  a 
foil.  These  were  unsuccessful,  as  the  material  work  hardened  to  such  an  extent  that 
it  could  not  be  flattened  to  any  appreciable  degree.  It  was  found  that  the  gold  could 
be  rendered  malleable  by  heating  at  150  °C  in  air  for  one  hour.  The  sample  was 
cut  into  several  smaller  pieces  with  a  ceramic  knife.  These  were  rolled  to  ~  0.1  mm 
thickness  using  the  mechanical  roller  in  the  Materials  Science  machine  shop.  In  order 
to  prevent  contamination  from  the  steel  rollers,  the  gold  was  sandwiched  between  two 
sheets  of  molybdenum4  of  three  9's  purity.  The  gold  sheets  obtained  from  this  rolling 
were  then  cut  into  16  foils  3  mm  wide  and  50  mm  long  with  a  pair  of  ceramic  scissors. 

The  heat  treatment  of  these  foils  was  straightforward,  as  it  was  found  that  RRR's 
of  over  1000  could  be  obtained  by  annealing  in  air  for  2.5  hours  at  700  °C.  The 
sample  was  annealed  for  3  days  at  900  °C  under  an  oxygen  pressure  of  10-4  torr. 
This  procedure  yielded  a  RRR  of  1700.  The  foils  were  welded  to  the  bottom  of  the 
silver  cross  with  an  acetylene  torch.  After  welding,  the  foils  were  dusted  with  silicon 
oxide  powder  to  prevent  them  from  sticking  together. 

B.3    Experimental  Search 

The  procedure  outlined  in  chapter  2  for  the  silver  experiment  was  used  in  the 
search  for  the  gold  NMR  signal.  The  absence  of  a  signal  led  to  a  very  large  search 
in  frequency,  however.  A  large  number  of  sweeps  were  made  at  frequencies  in  the 

3NIST,  US  Department  of  Commerce,  Gaithersburg,  MD  20899. 
4Schwarzkopff  Technologies  Corp.,  35  Jeffrey  Ave.,  Holliston,  MA  01746. 
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region  of  1  kHz  for  various  static  fields.  As  no  gold  signal  was  observed  in  this  region, 
several  attempts  were  made  to  observe  the  NMR  signal  by  demagnetizing  the  second 
magnet  straight  to  zero-field,  and  searching  for  a  resonance  from  10  to  100  Hz.  In  this 
regime,  no  signal  could  be  found.  For  the  second  gold  experiment  only,  an  attempt 
was  made  to  measure  the  equilibrium  magnetization  at  3  kHz  at  a  temperature  of 
100  /jK.  Unfortunately,  no  NMR  signal  was  detected  in  this  scheme,  either.  For  the 
time  being,  the  observation  of  the  C W  NMR  signal  in  gold  remains  an  open  topic. 
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